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Q ; ABSTRACT 

■ Aims. We have performed a comprehensive analysis of a sample of 20 starburst galaxies that show the presence of a 

^->' substantial population of very young massive stars, most of them classified as Wolf-Rayet galaxies. 

C^ ' Methods. In this paper, the last of the series, we analyze the global properties of our galaxy sample using multiwavelength 

^^ , data extracted from our own observations (Ha fluxes, B and //-band magnitudes) and from the literature, that include 

^^ . X-ray, FUV, FIR, and radio (both Hi spectral line and 1.4 GHz radio-continuum) measurements. 

lf> • Results. The agreement between our Ha-based star-formation rates (Sfr) and those provided by indicators at other 

wavelengths is remarkable, but we consider that the new Ha-based calibration provided by Calzetti et al. (2007) should 

be preferred over older calibrations. The FUV-hased Sfr provides a powerful tool to analyze the star-formation activity 

in both global and local scales independently to the Ha emission. We provide empirical relationships between the ionized 

gas mass, neutral gas mass, dust mass, stellar mass, and dynamical mass with the _B-luminosity. Although all mass 

estimations increase with increasing luminosity, we find important deviations to the general trend in some objects, that 

(~| , seem to be consequence of their particular evolutionary histories. The analysis of the mass-to-light ratios give similar 

O 11 results. We investigate the mass-metallicity relations and conclude that both the nature and the star-formation history 

I ' are needed to understand the relationships between both properties. The majority of the galaxies follow a Schmidt- 

O ' Kennicutt scaling law of star-formation that agrees with that reported in individual star-forming regions within M 51 

,y I but not with that found in normal spiral galaxies. Dwarf galaxies seem to be forming stars more efficiently than the 

C/2 , outskirts of spiral galaxies. We found a relation between the reddening coefScient and the warm dust mass indicating 

, ^/ that the extinction is mainly internal to the galaxies. The comparison with the closed-box model also indicates that 

environment effects play and important role in their evolution. 

I , Conclusions. Considering all multi- wavelength data, we found that 17 up to 20 galaxies are clearly interacting or 

K^ . merging with low-luminosity dwarf objects or H I clouds. The remaining three galaxies (Mkn 5, SBS 1054-1-364, and 

/^ ' SBS 1415-f437) show considerable divergences of some properties when comparing with similar objects. Many of the 

1^ ' interacting/merging features are only detected when deep optical spectroscopy and a detailed multi- wavelength analysis, 

vy—j I including H I observations, are obtained. We conclude that interactions do play a fundamental role in the triggering 

j.-^ , mechanism of the strong star-formation activity observed in dwarf starburst galaxies. 

I/-V , Key words, galaxies: starburst — galaxies: interactions — galaxies: dwarf — galaxies: abundances — galaxies: kinematics 

^— ^ ■ and dynamics — stars: Wolf-Rayet 

o: 

^ ' 1. Introduction ies have revealed the fossil remnants of interaction/merger 

• 1— I activity, increasing the evidences that interactions and 

*^ 1.1. Galaxy interactions and starburst activity mergers trigger star-formation phenomena in spiral 

^ . Since the discovery of the starburst galaxies Salaries (|Koribalski| MM iKennicut^ [Ml | Nikolic et al. | 

(iSargent fc Searld ITOTOl) . many studies have tried to ^^Q^)- Infrared observations confirmed the existence of 

understand the processes that trigger the strong star- ^^ry intense starbursts m major disk-disk mergers (e.g., 

formation activity in these objects. The hypothesis that J^^^P^ ^ Wright 1985; Solomon & Sage 1988; Sanders 

gravitational interaction (not necessary merging) of galax- ^ Mi'^abel 1996; Genzel et al. 1998; Arribas et al. 2004) 

ies enhances star formation or leads to starburst activity Actually, almost 100% of galaxies with far-mfrared {FIR) 

was made so on after the reco gnition of the starburst luminosities of aboujlO Lq are m interacting/merging 

phenomenon. iLarson fc Tinsl^ (IT978I) did a study of s^^tems ( |Sanders| ll997j. Furthermore, analysis of large 



o 
u 



normal and peculiar iAr^lmt) sample of galaxies and S^kxy surveys (e.g., CfA2: Barton Ge Her & Kenyon 2000; 
demonstrated that recent (< lO^ yr) star-formation is 2dF: Lambas et al. 2003; SDSS: Nikohc et al. 2004) has 
more likely to occur in interacting than in non-interacting Provided new evidences for interaction-induced starburst 



galaxies. Since then, numerous studies of individual galax- 



activity. 



According to hierarchical clustering models of galaxy 
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1991; Kauffman & White 1993; Springer et al. 2005). 
Observations of local and distant luminous blue galaxies 
(LBG) and Lyman break galaxies seem to confirm that 
galaxy interactions are more common at high redshifts. 
(e.g., Guzman et al. 1997; Hopkins et al. 2002; Erb et 
al. 2003; Werk, Jangren & Salzer 2004; Colina, Arribas 
& Monreal-Ibero 2005; Overzier et al. 2009; Cardamone 
et al. 2009) but many details are still unclear (i.e., Basu- 
Zych et al. 2009). Indeed, detailed studies of local interact- 
ing/merging galaxies provide vital clues to our knowledge 
in galaxy formation and evolution, as they constrain the 
properties of the hierarchical formation models. 

Recent observations also suggest that interactions and 
mergers between dwarf galaxies also trigger the star- 
formation activity and play a fundamental role in the evo- 
lution of dwarf galaxies (i.e., Mendez & Esteban, 2000; 
Ostlin et al. 2001, 2004; Bergvall & Ostlin 2002; Johnson 
et al. 2004; Bravo-Alfaro et al. 2004, 2006; Gumming et 
al. 2008, Garcia-Lorenzo et al. 2008; Lopez-Sanchez &; 
Esteban, 2008, 2009; James et al. 2010). Many of these 
studies have been done in Blue Compact Dwarf Galaxies 
(BCDGs), that are low-luminosity, low-metallicity (^10 % 
solar) galaxies showing compact and irregular morphologies 
and undergoing an intense and short-lived episode of star 
formation (i.e., Izotov & Thuan 1999; Cairos et al. 2001a, b; 
Papaderos et al. 2006), on top of an old underlying popu- 
lation with ages of several Gyrs (Noeske et al. 2003, 2005; 
Amo rin et al. 2007, 2009). Recent numerical simulations 
(jBekki 2008) satisfactory explain the physical properties 
of BCDGs as a consequence of the merging of two dwarf 
galaxies with larger fraction of gas and extended gas disks. 

Actually, much of our knowledge in interacting galaxies 
has been provided via H i observations. Neutral hydrogen 
gas is the best tracer for galaxy-galaxy interactions because, 
such the H i distribution is usually several times larger than 
the optical extent, it is more easily disrupted by exter- 
nal forces (tidal interactions, gas infall, ram pressure strip- 
ping) than the stellar dis k (jBroeils fc van WoerdenI [19941 : 
ISalpeter fc Hoffman! 1 1996D . The distribution and kinemat- 
ics of atomic gas within galaxies usually is more or less 
regular, but in many cases they revealed complex entities 
between galaxies such as tails, ripples and bridges, arcs, 
or independent Hi clumps that, in many cases, show lit- 
tle disturbance in their corresponding optical images (e.g., 
Schneider et al. 1989; Yun, Ho & Lo 1994; Hibbard & van 
Gorkom 1996; Verdes-Montenegro et al. 2001, 2002, 2005; 
Putman et al. 2003; Koribalski et al. 2003; 2004; 2005; 
Temporin et al. 2003, 2005; Emonts et al. 2006; Ekta et 
al. 2008; Koribalski & Lopez-Sanchez 2009; English et al. 
2010; see also The Hi Rogues Gallery, Hibbard et al. 2001). 
Several interferometric H i surveys, such as The H i Nearby 
Galaxy Survey (THINGS, Walter et al. 2008), the Local 
Volume H i Survey (LVHIS, Koribalski 2008) or the Faint 
Irregular Galaxies GMRT Survey (FIGGS, Begum et al. 
2008), are nowadays providing accurate H i and dynamical 
masses in hundreds of nearby galaxies, many of them be- 
ing dwarf objects, as they acc ount for ^^85% of the kno wn 
galaxies in the Local Volume (jKarachentsev et al.ll2004l) . 

1.2. The interplay between gas and stars in galaxies 

However, to understand interaction processes in dwarf 
galaxies we first have to know how stars and gas inter- 



act in low-mass environments. Indeed, feedback from mas- 
sive stars is the dominant process that affects the interstel- 
lar medium (ISM) of these galaxies. Violent star-formation 
phenomena may disrupt the galaxy's gas and even expel 
it to th e intergalactic medium, as some theoretical models 
predict ()Mac Low fc Ferraralll999() . But alternative models 
(e.g., Silic h fc Tenorio-T agle 1998) and the available ob- 
servations (lBoman£ll2005[ ) suggest that dwarf galaxies keep 
their processed material. Furthermore, the links between 
the observational characteristics (fluxes, colors, morpholo- 
gies or sizes) and the underlying physical properties of the 
galaxies (stellar, dust, gas, baryonic, and dark matter con- 
tent, chemical abundances, star formation rate, star forma- 
tion history) are still not well known. 

For example, there are still many caveats in the under- 
standing of the interplay between the star formation rate 
(Sfr) and the properties of the ISM. A very important 
step was achieved with the Schmidth-Kennicut t power-law 
relation (|Schmidtlll959l 119631 : lKennicuttlll998D that corre- 
lates the average Sfr per unit area and the mean surface 
density of the cold gas (atomic plus molecular). But trac- 
ers of star-formation, including optical colors and Ha flux 
(e.g., Larson fc Tir isley, 1978: Kennicutt 1998; Ca lzetti et 
al. 2007), FIR fl ux (lKemiicutt|l998t lHeckmanlll999[ ). radio- 
continuum flux ('Condon"1992^, and far-ultraviolet (FUV) 
flux (Kennicutt 1998; Salim et al. 2007), often yield to very 
different values of the Sfr. Although the density of atomic 
gas is known in some cases, no many direct measurements 
of the molecular gas are available, being specially rare in 
dwarf galaxies (i.e., Taylor et al. 1998; Barone et al. 2000; 
Braine et al. 2000, 2001, 2004). 

On the other hand, the physics underlying the relation- 
ship between stellar mass (or luminosity) with the metal- 
licity is still far from clear, besides the important observa- 
tional (e.g., Tremonti et al. 2004; van Zee & Haynes, 2006; 
Kewley fc Ellison 2008) and theoretical (e.g., De Lucia et 
al. 2004; Tissera et al. 2005; De Rossi et al. 2006; Dave fc 
Oppenheimer 2007) efforts that aimed to explain it. Indeed, 
one of the main problems is to derive the real metallicity 
of the ionized gas, as empirical calibrations based on the 
direct estimation of the electron temperature (Te) of the 
ionized gas and theoretical methods based on photoioniza- 
tion models provide very different oxygen abundances (e.g., 
Yin et al. 2007; Kewley fc Ellison 2008; Esteban et al. 2009; 
Lopez-Sanchez fc Esteban, 2010). Finally, the present un- 
derstanding of correlations between the H i content, stellar 
populations and star formation in dwarf starburst galaxies 
is still at a preliminary stage because of the lack of detailed 
optical/ A^/i? images and spectra and/or interferometric H i 
maps in these systems. 



1.3. About this paper series 

In our paper series, we have presented a detailed pho- 
tometric and spectroscopic study of a sample of strong 
star-forming galaxies, many of them previously classified 
as dwarf galaxies. The majority of these objects are Wolf- 
Rayet (WR) galaxies, that are a very inhomogeneous class 
of star-forming objects which have in common that the 
ongoing or most recent star formation event has pro- 
duced stars sufficiently massive to evolve to the WR stage 
(Schacrcr, C ontini fc Pindao 1999 ). The presence of WR 
features in the spectra of a galaxy constrains the properties 
of the star-formation processes. As the first WR stars typ- 
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ically appear around 2-3 Myr after th e starburst is initi- 
ated and disappear within some 5 Myr (jMevnet &: Maeden 
l2005f) . their detection informs about both the youth and 
strength of the burst, offering the opportunity to study 
an approximately coeva l sample of very young starbursts 
(jSchaerer fc Vaccalll998D . 

Our main aim is the study of the formation of mas- 
sive stars in starburst galaxies and the role that interaction 
with or between dwarf galaxies and/or low surface bright- 
ness objects have in its triggerin g mechanism. In Paper I 
(jLopez-Sanchez fc Estebanll2008D we exposed the motiva- 
tion of this work, compiled the list of the analyzed WR 
galaxies (Table 1 of Paper I) and presented the results of 
optical/JV/fi broad-band and Ha photometry. In Paper II 
(jLopez-Sanchez fc Estebanll2009f ) we presented the results 
of the analysis of intermediate-resolution long slit spec- 
troscopy of 16 objects of our sample of WR galaxies -the 
results for the other 4 objects have been pubHshed sep- 
arately. In Paper II, we also specified the oxygen abun- 
dances of the ionized gas (they were computed following 
the direct Tg method in the majority of the cases ) and 
analyzed the kinematics of the io nized gas. In Paper III 
(|L6Dez-Sanchez fc Estebanll2010a[ ). we studied the O and 
WR stellar populations within these galaxies, and com- 
pared with theoretic al evolutionary synthesis models. In 
Paper IV (Lope z-Sanchez fc Estebanl|2010b ), we analyzed 
globally the optical/iV7i? properties of the galaxies, con- 
cluding that such detailed analyses are fundamental to un- 
derstand the star-formation histories of the galaxies. In this 
paper, the last of the series, we perform a comprehensive 
multiwavelength analysis considering all the optical and 
NIR data but also including radio, FIR, FUV and X-ray 
data available in the literature. 

The selection criteria of the galaxy sample were the fol- 
lowing. We used the most rec ent catalogue of Wolf-Rayet 
galaxies (jSchaerer et al.lll999f ). which contains a very in- 
homogeneous group of starbursting objects, to make a 
list of dwarf objects that could be observed from the 
Northern Hemisphere. Hence, we did not consider either 
spirals galaxies or giant H ii regions within them, and con- 
sidered only dwarf objects, such as apparently isolated 
BCDGs and dwarf irregular galaxies that had peculiar mor- 
phologies in previous, shallower imaging. We also chose 
two galaxies belonging to the Schaerer et al. (1999) cat- 
alogue that were classified as suspected WR galaxies (Mkn 
1087 and Tol 9), to confirm the presence of massive stars 
within them (see Papers II and HI). The galaxy IRAS 
08339-1-6517 was also included because previous multiwave- 
length results suggested that the WR stars could still be 
present in its youngest star-forming bursts (see Lopez- 
Sanchez et al. 2006). With this, we got a list of '^40 sys- 
tems to observe and analyze using the telescopes avail- 
able at Roque de los Muchachos (La Palma, Spain) and 
Calar Alto (Almeria, Spain) observatories. We added the 
southern galaxy NGC 5253, for which we obtained deep 
echelle spectrophotometry using 8.2m VLT, because of the 
very intriguing properties it possesses (see Lopez-Sanchez 
et al. 2007, 2010). The final sample of 20 galaxies was cre- 
ated considering those galaxies for which we obtained opti- 
cal/ Affl broad-band and Ha images plus the deep optical 
spectroscopy during our observation runs. We already have 
all these data for other '^15 galaxies, the analysis of these 
systems will be presented in the future elsewhere, but its 
preliminary results seem to agree with the main results re- 



ported in this paper. Hence, our galaxy sample is not com- 
plete, but we consider it represents quite well dwarf galax- 
ies experiencing a very strong star- formation burst. Indeed, 
this was the main bias introduced when choosing the galaxy 
sample, such as we focused only in galaxies in which WR 
stars are detected. It would be very interesting to extend 
this analysis to similar star-forming galaxies that do not 
show WR features, as the samp l e of B CDGs analyzed by 
iGil de Paz. Madore fc Pevunoval (|2003[ ). 

The structure of this paper is the following. In Sect. 2 
we describe the details of the radio, FIR, FUV, and X-ray 
data extracted from the literature, providing some very use- 
ful relations. Sect. 3 analyzes the star-formation activity in 
our sample galaxies considering all multi-wavelength cali- 
brators to the Sfr. We check if our sample galaxies follow 
the radio/ FIR correlation in Sect. 4. Next, Sect. 5 com- 
piles, analyses and compares all mass estimations derived 
in this work. Several mass-metallicity relations are inves- 
tigated in Sect. 6. We study whether our galaxies satisfy 
the Schmidt-Kennicutt relation in Sect. 7. Section 8 ana- 
lyzes and compares several mass-to-light ratios. The dust 
properties within our starburst galaxies are investigated in 
Sect. 9. We compare the predictions of the closed-box model 
with our observational data in Sect. 10. Finally, Sect 11 
compiles a quantitative analysis of the interaction features 
considering all available multi-wavelength data. The con- 
clusions reached in our analysis are compiled in Sect. 12. 
The Appendix describes the main results found in each of 
the analyzed WR galaxies. 

Hence, this is essentially an observational work. Each 
system has been carefully analyzed considering all avail- 
able data (those specifically obtained for this work and 
those compiled from literature) with the final aim to un- 
derstand its chemical and dynamical evolution, its stellar, 
dust, gas, and dark matter content, the relative importance 
of its stellar populations (WR, young, intermediate- age and 
old stars) and its star formation properties. Our data sup- 
port the hypothesis that interactions between galaxies and 
dwarf or low surface-brightness objects (that can not be de- 
tected using less detailed and less deeper observations) have 
a considerable importance in the triggering mechanism of 
massive star formation activity in this kind and young star- 
bursts. We have produced the most complete, detailed, and 
exhaustive data set of this kind of galaxies, so far, involving 
multi-wavelength data and a careful analysis of each indi- 
vidual object following the same procedures and equations. 



2. Multi-wavelength data completeness 

We have performed an exhaustive literature search to com- 
plete the optical/ A'/i? observations of our WR galaxy sam- 
ple with data from other wavelengths (radio, far-infrared, 
far-ultraviolet, and X-ray). Here we describe all these data 
and the useful properties we derive from them. 



2.1. Radio data 

2.1.1. HI data at 21 cm 

Observations in the hyperfine transition of the neutral hy- 
drogen, H I, with a rest frequency of 1420.405 MHz, have 
been key to understand the distribution and kinematics of 
the atomic gas within galaxies, included the Milky Way. 
Neutral gas observations are very important because they 
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Table 1. Radio data compiled from the literature for our WR galaxy sample. We include the flux of the 21 cm H I emission line, 
-Fhi, its equivalent width, Wui, and the radio-continuum flux at 1.4 GHz, S'i.4ghz. 



Galaxy 



[Jy km s^-"-] 



[km s^-*-] 



Ref. 



Si. 4 GHz 

[mJy] 



Ref. 



HCG 31 AC 


5.15 


169.2A+190.6C 


VM05 


22+3 


VM05 


HOG 31 B 


2.74 


85.8 


VM05 


2.1+0.3 


VM05 


HCG 31 F 


0.866 


74.6 


VM05 






HCG 31 G 


2.74 


84.9 


VM05 


3.3+0.5 


VM05 


Mkn 1087 


5.38 


270 


GG81 


12.1+0.6 


Co98 


Haro 15 


3.11+1.01 


220 


GG81 


17.8+1.0 


Co98 


Mkn 1199 


1.78+0.67 


170 


DC04 


36.2+1.2 


Co98 


Mkn 5 


2.12+0.27 


22.4+4.9 


Pa03 


<2.8 


HSLD02 


IRAS 08208+2816 








15.2+0.6 


Co98 


IRAS 08339+6517 


3.68+0.46 


~300 


Ca04 


33.56"^ 


Co90 


POX 4'' 


4.31 


130 


LSlOb 


4.2+0.5 


Co98 


UM 420 








1.1+0.3 


HSLD02 


SBS 0926+606A<= 


1.30+0.49 


120+37 


P02 


2.7+0.6 


HSLD02 


SBS 0926+606B^ 


1.10+0.49 


120+37 


P02 






SBS 0948+532 








<0.9 


HSLD02 


SBS 1054+365 


4.03+0.39 


117+11 


ZOO 


1.28+0.14 


BWH95 


SBS 1211+540 


0.71+0.12 


47 


H05 


<0.9 


HSLD02 


SBS 1319+579 


8.4 


134 


H07 


<2.9 


HSLD02 


SBS 1415+437 


4.73+0.32 


66 


H05 


<0.5 


H05 


HI Zw 107 


4.48+0.79 


200+25 


P03 


8.0+0.5 


Co98+Y01 


Tol9^ 


5.02+0.40 


185 


LSlOb 


19.2+0.7 


Co98 


Tol 1457-262'* 


4.3 


176 


Kor06 


38.9+1.8 


Co98+Y01 


Arp 252'''<= 






Kor06 


97.6+3.0 


Co98+Y01 


NGC 5253^ 


43.1+2.6 


106+6 


LSlOa 


87.1+3.5 


LSlOa 



"^ Co90 llCondon et al.lll99CD gave radio-continuum valu es at 1.49 GHz. The value of Si.4GHz shown in the table was computed from 
S'l 49 GHz using the ("Condon. Cotton fc Bro derickll2002D relation between both quantities, as it was explained in 

(200l). 



I Lope z-Sanchez. Esteban & Garcia-Rojas 

' A detailed analysis of the H I gas within the galaxies POX 4, Tol 9 and NGC 5253 using the Australia Telescope Compact Array (ATCA) 
will be soon presented elsewhere lUjopez-Sanchez ct al. 2010a b). 

'^ SBS 0926+606 was observed bv iHuchtmcier et al.i (i2007i ) . who gave a measurement of the Hi flux for both A and B galaxies. Only 
interferometric studies can disentangle the a mount of neutral gas i n the individual galaxies. 

■* This galaxy was observed in Hi bv iCasasoJa. Bettoni fc Galletal | |2004D using a single-dish antenna, but it was not detected. 
'^ This galaxy is not detected in Hi in HIPASS (Koribalski 2006, priv. comm.). 

Cannon et al.l 



Galletal | [200l ): Co90:rCon don ct al.' (1990): Co98: 'Condon et al.l 
GG81: lGordon fc G ottcsman ( 198ll): HSLD02: iHopkiiis. Schultc-Ladbcck & Drozdovskv (2 003) : 



Refer ences: Ca04: 'Can non et al.l ll2004f) : CBG0 4: [ Casasola. B ettoni 

(|1998|1: DC04: f Pavoust fc Continil ll2004fl: GG 81: lGordon fc G ottcsm _ 

H05: IHuchtmcier. Krishna & Pctrosian (2005") ;H07: Huchtmeicr et al. (2007); Kor06: Koribalski (2006), priv. com m.; LSlOa: 
ILopez-San chez et al. (2010a); LSlOb: Lopez-Sanchez ct al. (2010b); Pa03: Paturcl et al. ( 20 0j) : Hype rLEDA; P02: IPustiInik et al.l (|2002l ): 
VM05: IVerdes-Montenegro et al.1 II2005I) : YOl: lYun. Reddv fc CondonI <2001l) : ZOO: IZasov et ahl (l2000fl . 



are used to determine both the neutral gas mass (H i gas) 
and the dynamical mass (Mdyn) of the systems. Single-dish 
H I surveys, -e.g. Mathewson et al. 1992, the H i Parkes Sky 
Survey (HIPASS, Barnes et al. 2001; Koribalski et al. 2004; 
Meyer et al. 2004), and the Arecibo Legacy Fast ALFA sur- 
vey (ALFALFA; Giovanelli et al. 2005)-, give spectra with 
detected H i emission of thousands of galaxies. However, the 
best tool to analyze the neutral gas content in galaxies is via 
radio interferometer observations (e.g., THINGS; LVHIS; 
FIGGS; The Hi Rogues Gallery). Knowing the amount of 
available neutral gas, the timescale of the starbursts (i.e., 
the time in which the Hi cloud will be exhausted if the 
star formation activity continue at the current Sfr) can be 
calculated. 

Table [1] compiles all H i 21 cm data found for our galaxy 
sample. The majority of the H i data is provided by single- 
dish Hi observations, but for some few cases (HGC 31 and 
IRAS 08339+6517) interferometric H i maps are available. 
Table [1] lists the H i flux density, /h i (in units of Jy km 
s""'^), and the Hi equivalent width, Whi (in kms^^). We 
note that for 3 galaxies (POX 4, Tol 9 and NGG 5253) 
we are using the data provided by our new interferomet- 
ric maps obtained using the Australia Telescope Compact 



Array. For these objects, we compile the integrated H I 
flux and width; their detailed analysis will be soon pre- 
sented elsewhere (JLopez-Sanchez et al.ll2010al lbr). The total 
H I mass is computed applying 

Mhi = 2.356x lO^d^Hi (1) 

(|Robertdll975HRoberts fc Havneslll994{ ) where the distance 
to the galaxy, d, is expressed in Mpc and the result for the 
neutral gas mass is given in solar units. The dynamical mass 
of the system, Mdyn, can be estimated from H i radio ob- 
servations considering the inclination-corrected maximum 
rotation velocity, WmaKi that is obtained at radius Rmax 
and assuming a virial equilibrium. 



Mdy^ = 2.31 X IQ'^R^axivl^ax?. 



(2) 



Wui 



being the result in solar masses and assuming f L„.^ — „ . . . 

. . . . ^ max 2 sin i 

The inclination angle, i, is defined as that found between 
the plane of the sky and the plane of the galaxy (hence, 
1=90° in an edge-on galaxy and i—0° in a face-on galaxy). 
We usually estimated this angle assuming that the elliptical 
shape of the galaxy is just a consequence of its orientation. 
Note that the usual problem deriving virial masses is the 
unknowledge of the inclination angle, i, and sometimes also 
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Table 2. FIR and FUV data for the WR galaxy sample analyzed in this work. FIR data were provided by IRAS, FUV data 
were provided by GALEX. The value of the total FIR flux was computed using Equation [7] 



Galaxy 


/l2 (jm 


/25 Mm 


/eo ^lm 


/lOO fim 


Ffir 


Reg. 


mpuv 


Ifuv 




[Jy] 


[Jy] 


[Jy] 


[Jy] 


(a) 


(b) 


[mag] 


(c) 


HOG 31 


0.110+0.020 


0.580+0.040 


3.92+0.31 


5.84+0.47 


2.01+0.16 


AC 
B 
E 
F 
G 
H 


14.97+0.06 
16.80+0.09 
18.01+0.08 
18.24+0.11 
15.93+0.06 
20.63+0.40 


7.54+0.44 

3.32+0.26 

0.456+0.035 

0.460+0.047 

3.11+0.18 

0.0041+0.0014 


Mkn 1087 


0.103+0.029 


0.414+0.058 


3.03+0.33 


4.44+0.40 


1.54+0.16 








Haro 15 


0.118+0.034 


0.297+0.089 


1.36+0.12 


1.97+0.20 


0.690+0.064 




15.12+0.07 


9.41+0.57 


Mkn 1199 


0.282+0.031 


1.28+0.09 


6.82+0.34 


8.85+0.53 


3.33+0.18 








Mkn 5 


<0.0503 


<0.0533 


0.21+0.04 


<0.8473 


<1.75 




17.76+0.07 


1.59+0.10 


IRAS 08208+2816 


0.126+0.029 


0.278+0.067 


1.15+0.09 


1.70+0.17 


0.588+0.051 




16.59+0.07 


3.75+0.23 


IRAS 08339+6517 


0.250+0.025 


1.13+0.02 


5.81+0.04 


6.48+0.09 


2.71+0.02 


c 


15.37+0.06 
19.70+0.10 


10.7+0.6 
0.160+0.014 


POX 4 


<0.987 


0.153+0.040 


0.629+0.057 


<0.5798 


<0.278 


POX 4 
Comp. 


16.03+0.06 
19.06+0.12 


2.82+0.16 
0.269+0.030 


UM 420 




<0.275 


0.411: 


0.613: 


0.211: 




18.22+0.06 


0.377+0.021 


SBS 0926+606 


<0.07553 


<0.08818 


0.269+0.046 


<0.5296 


<0.154 


A 
B 


16.71+0.06 
17.84+0.10 


1.75+0.09 
0.829+0.073 


SBS 0948+532 


















SBS 1054+365 


<0.055 


<0.100 


0.536+0.048 


0.97+0.15 


0.296+0.035 




16.67+0.07 


1.18+0.08 


SBS 1211+540 














18.40+0.07 


0.370+0.022 


SBS 1319+579 






0.209: 


0.685: 


0.154: 


A 


17.00+0.06 
18.60+0.06 


0.864+0.048 
0.199+0.010 


SBS 1415+437 














16.22+0.06 


3.95+0.22 


III Zw 107 


<0.0968 


0.336+0.050 


1.37+0.20 


1.72+0.31 


0.662+0.104 








Tol 9 


0.111+0.030 


0.465+0.051 


2.71+0.22 


<5.516 


<1.58 




17.57+0.07 


4.22+0.22 


Tol 1457-262 


<0.117 


0.611+0.067 


3.09+0.19 


3.68+0.40 


1.47+0.11 


Obj 1 
Obj 2 

#15 


15.84+0.07 
16.15+0.06 
17.41+0.13 
19.71+0.21 


6.95+0.47 

5.26+0.29 

1.64+0.18 

0.199+0.037 


Arp 252 


0.188+0.023 


0.994+0.050 


3.91+0.20 


4.11+0.25 


1.79+0.10 


A 
B 


18.67+0.08 
18.84+0.07 


1.91+0.13 
0.441+0.027 


NGC 5253 


2.50+0.02 


12.07+0.05 


29.84+0.07 


30.08+0.21 


13.49+0.05 




12.81+0.06 


123.0+6.5 



(a) In units of 10 ^^ erg s ^ cm ^. 

(b) Region within each system (HCG 31, IRAS 08339+6517, POX 4, SBS 0926+606, SBS 1319+579, Tol 1457-262, and Arp 252; see Paper I 
for identification of the regions). The FIR emission provided by IRAS does not allow to distinguish between these regions, but FUV data 
provided by GALEX does. In Arp 252, region A is galaxy ESO 566-8 and region B is galaxy ESO 566-7. 

(c) In units of 10" ""^^ erg s~^ cm~^ A"-*^. 



Umax, specially in galaxies showing disturbed morpholo- 
gies. We adopted the maximum radius observed in our deep 
optical images. Therefore, as the extension of the neutral 
gas is usually larger than the extension of the stellar com- 
ponent, our values of Mdyn may b e underestimated. The 
gas depletion timescale defined by ISkillman et al.l (|2003[ ) 
was computed using A/hi and the assumed Sfr derived for 
each galaxy (see below). 



2.1.2. Radio-continuum data 

For an individual star-forming galaxy, the Sfr is directly 
proportional to its radio luminosity (i.e., Condon 1992). 
Hence, the radio continuum flux is widely used as a dust- 
free indicator of the star formation rate. Nearly all of the 
radio-continuum luminosity from galaxies without a sig- 
nificant Active Galactic Nucleus (AGN) can be traced to 
recently formed massive (M > 8 Mq) stars (Condon et 
al. 1992). The 10% of the continuum emission at 1.4 GHz 
is due to free-free emission from extremely massive main- 
sequence stars (thermal emission) and almost 90% is syn- 
chrotron radiation from relativistic electrons accelerated 
in the remnants of core-collapse supernovae (non-thermal 



emission). As the stars that contribute significantly to the 
radio emission have lifetimes r < 3 x 10^ yr and the 
relativistic electrons have lifetimes r < 10* yr, the cur- 
rent radio luminosity is nearly proportional to the rate 
of massive star formation during t he past r < 10* yr 
(|Condon. Cotton fc Broderickll2002[ ): 



SFRiAGU. [M > 5Mq) ^ 2.5 x 10-""Li.4ghz 



(3) 



where Li.4ghz has units of W Hz~^. 

Table [T] compiles all 1.4 GHz radio-continuum flux data 
available for our WR galaxy sample in the literature. The 
1.4 GHz luminosity . iS'i.4ghzi can be comp i ited u sing the 
expression given bv lYun. Reddv fc CondonI ()2001[ ): 



log Li. 4 GHz == 20. 07 -I- 2 log d-t- logs'!. 4 GHz, 



(4) 



where the result is given in units of W Hz"^, the distance 
d is expressed in Mpc and 151.4 ghz is expressed in Jy. 

Radio-continuum observations at several cm wave- 
lengths are used to quantify the thermal and non-thermal 
contributions, and thereby distinguish older and supernova- 
rich regions from younger and mostly thermal areas 
(i.e., Deeg et al. 1993, Beck et al. 2000, Cannon et 
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al. 2004, 2005). These observations also permit to de- 
tect extremely young, dense heavily embedded sta r clus- 
ters pCo bulnickv & Johnson 1999; Johnson & Kob ulnickvl 
l2003f) . Although radio data at frequencies different of 
1.4 GHz are not usually available for this kind of galaxies, 
we applied the equation provided bv lDopita et al.l ([2002:), 



^1.4 GHz thermal — 1-21 X 10 Fffa, 



(5) 



to obtain an estimation of the thermal emission at 1.4 GHz, 
-Pi.4 GHz thermal, usiug the Ha flux derived from our im- 
ages (see Paper I). In this equation, Fho is in units of 
erg cm~^ s~^ and the result is given in mJy. The com- 
parison between -Fi.4ghz thermal and Fi.4ghz allows the 
estimation of the non-therm al flux. Gondon (1992) and 
iNiklas. Klein fc Wielebinskil (|T997) indicated that the non- 
thermal component is more than 90% of the total at this fre- 
quency. It is common to conside r the non-thermal to ther- 
mal ratio, i?: lDopita et al.l (|2002h reported that the average 
value in starburst galaxies is logi? = 1.3 ± 0.4. Radio con- 
tinuum and FIR data help to discern between the normal 
or active (that is, a galaxy hosting an AGN) nature. 

2.2. FIR data 

Many of the problems found to derive the Sfr from optical 
data can be avoided by measuring the far-infrared (FIR) 
and sub-millimeter spectral energy distributions (SEDs). 
These are determined by the re-radiation as thermal con- 
tinuum by the dust grains of stellar photospheric radiation 
absorbed in the visible and UV regions of the spectrum. 
Assuming that the dust completely surrounds the star form- 
ing regions, it acts as a bolometer reprocessing the luminos- 
ity produced by the stars. Therefore, the Sfr can be also 
computed using the oretical ste l lar flu x distributions and 
evolutionary models. [Kennicutg (|1998[ ) provides the follow- 
ing correlation between the Sfr (in units of Mq yr~^) and 
the far-infrared flux: 



SFRfib. = 4.5 X lO^^^'L 



FIR, 



(6) 



where Lpjji (given in units of erg s ^ ) i s obtained using the 
FIR&ux between 42.5 and 122.5 /Ltm (jSanders fc Mirabell 



Ffir = 1.26 X 10' 



(2.58/60 + /loo), 



(7) 



being /go and /loo the flux densities (in Jy) for 60/Ltm and 
100/Ltm and the conventional expression between flux and 
luminosity, L — AncPF. This relation can be applied only 
in starbursts with a ges less than 10^ y r, where the approxi- 
mations assumed bv lKennicuttI ()1998f ) are valid. If the SFR 
value derived from Lpjfi agrees with that estimated from 
the Ha luminosity we may consider that the correction by 
extinction done to derive the Ha flux is correct. 

Assuming that all the UV and blue radiation from mas- 
sive stars is absorbed by grains and is re-emitted as thermal 
radiation in the 40-120 ^m band, Gondon (1992) derives the 
following relation between Sfr and Lgo^m (in units of W 
Hz-i): 



(8) 



SFRgo^m '' 


- 1.96 X 10"^ 


-4 T 


i.e.. 






SFRgQ^m '' 


- 2.346 X 10" 


feO fj.md , 



being /eo um expressed in Jy a nd d in Mpc. Although some 
authors (jLonsdale et al.l [19841 ) have argued that radiation 
from stars less massive than 5 M© will contribute signif- 
icantly to FIR emission from galactic disks, this relation 
see ms to give good valu es for the SFR. 

iRoussel et al.l ()2001h provide an alternative Sfr calibra- 
tion using the 15 /im luminosity: 



SFR 



15 fl7n 



3.66 X 10-^r/i5 



fim ; 



(10) 



being fiSf^m the monochromatic flux at 15 /im and d the 
distance in Mpc. This formula is applicable only when the 
mid-infrared emission is dominated by unidentified infrared 
bands (UIBs) with a negligible very small grains (VSG) 
continuum, which is the case in disk galaxies, but is not 
always verified at 15 /zm in galactic central regions. It may 
be assumed that Li^^m ~ Li2^rn- 

The warm dust mass can be estimated using the 60 
and 100 iiva. fluxes and applying the relation given by 
iHuchtmeier. Sage fc Henke l (1995), 



Mdust — 4.78d /loo/im 



exp 



2.94 



V /fin ;*m / -I / 



60 /j-in 



(9) 



where the distance is expressed in Mpc, the flux densities 
are in Jy and the result is given in Mq. 

We have used the far-infrared (FIR) data provided by 
the Infrared Astronomical Satellite (IRAS) to obtain the 
monochromatic fluxes at 12, 25, 60 and 100 /im. These data 
are used to get an independent estimation of the Sfr and 
to derive the warm dust mass within every galaxy. We will 
also check if the galaxies follow the F/i?-radio relationship. 
Table [2] compiles all the FIR data found for our sample of 
WR galaxies, three of them have no useful measurements 
at these frequencies. 

2.3. FUV data 

In the last years, the GALaxy Evolution eXplorer (GALEX) 
satellite is providing astonishing ultraviolet {UV) images 
of galaxies, and revealing recent star- formation activity 
in their external regions (i.e., Gil de Paz, 2005, 2007; 
Thilker et al. 2005; Koribalski & Lopez-Sanchez, 2009). The 
GALEX point spread function in the central 0.5° has a full 
width at half-maximum [FWHM) of ^^5 arcsec, matching 
quite nicely with the spatial resolution of our optical/ A'7i? 
images. We searched for GALEX observations of the galax- 
ies that compose our sample in the far — C/V^-band {FUV, 
1350-1750 A), aU of them except four objects (Mkn 1087, 
Mkn 1199, SBS 0948-f 532 and III Zw 107) have useful FUV 
data. 

In general, the FUV emission of our sample galaxies 
matches quite well with their optical emission. In many 
cases, FUV emission is much more extended that the Ha 
emission. Figure [1] shows as examples the GALEX FUV 
images of HCG 31, Haro 15, and SBS 0926-f 606. As we can 
see when comparing with our optical images (see Paper I), 
the star-forming regions are clearly observed in FUV. Just 
quick comments about three galaxies: 

— we note that the eastern tail of SBS 0926+606 B is 
quite bright in the FUV image, suggesting a extended 
distributions of massive OB stars that we do not detect 
in our deep Ha image (see Fig. 15 and Sec. 3.10.1 in 
Paper I), 
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HCG31 



Q •% 






m 



A" 



GALEX FUV 



2' = 30.8 kpc 



HARD 15 






GALEX FUV 




2' = 100,8 kpc 



SBS 0926-1-606 






GALEX FUV 



2' = 32.4 kpc 



Fig. 1. Example of GALEX images, showing the FUV emission in HCG 31, Haro 15 and SBS 0926+606. Regions within each 
object have been labeled following the notation given in Paper I. 



— the star-forming galaxy #15 in Tol 1457-262 is clearly 
detected in the FUV emission, but the faint galaxy #16 
is not seen (Fig 31 and Sect. 3.18.1 in Paper I), and 

— in Arp 252, FUV emission is not only detected at the 
center of the galaxies (ESO 566-8 and ESO 566-7) but 
also throughout the tails and in tidal dwarf candidates 
c, e and d (see Fig. 34 and Sect. 3.19.1 in Paper I). 

FUV data can be used to get an independent estima- 
tion of the Sfr of the galaxies, f/l^-emission probes star 
formation over time-scales of ~100 Myr, the lifetime of 
the massive OB stars. We integrated the counts per sec- 
ond (CPS) within each galaxy or region a nd then app lied 
mpuv = -2.51og(CPS) -I- 18.82 (Morrissc v et all [20051) to 
derive the magnitude in the FUV-hand. We compile the 
rriFuv found for each galaxy in Table [2j We then corrected 
by extinction using Apuv = 7.9E{B — V) -the value of 
E(B — V) adopted for every region was derived from our op- 
tical spectroscopy and it is compiled in Table 5 of Paper I- 
and then apply 



Table 3. X-ray data available for our WR galaxy sample. 



fpuv = 1-40 X IQ-^'" X 10' 



-,0.4(18. 82-mO.^^,) 



(12) 



where rn'puy = mpuv — Apuv and fpuv is obtained in 
units of erg s^^ cm^^ A~^. The value of fpuv computed 
for each galaxy is also shown in Table [21 Once the FUV 
luminosity is computed (Lpuv = ^T^d^ Ifuv), the FUV- 
based Sfr i s deriv ed applying the calibration provided by 
iSalim et all (l2007h . 



SFRpiJY — 8.1 X 10 Lpjjy 



2.4. X-ray data 



(13) 



Finally, we also looked for the X-ray data available for 
our WR galaxy sample. Only four objects (HCG 31 AC, 
IRAS 08339-^6517, Tol 9 and NGC 5253) have been ob- 
served at these high frequencies; their X — ray luminosi- 
ties are compiled in Table [31 Beside these data, we will 
also use (see Section 3.3) the WR galaxy sample that 
[Stevens &: Strickland! (|1998a ,b) observed in X-ray. 

3. Analysis of the star formation rates 

The star formation rate (Sfr), defined as the stellar mass 
formed per unit time, is the standard parameter used to 



Galaxy 



logL; 



X,(0.2-2.0kcV) 



[ergs 



Rcf. 



HCG 31 


40.88 ± 0.13 


IRAS 08339-1-6517 


41.45 


Tol 9 


<40.29'' 


NGC 5253 


38.60 ± 0.18 



SS98 

SS98 

FFZ82 

SS98 



" This value was derived assuming the distance to Tol 9 and the flux 
quoted by FFZ82, /o.5-3kcV < 1-2 x lO'^^ ergcm-2s-\ and 
multiplying for 0.72 to correct for the X-ray range. 
References: FFZ82: Fabbiano, Feigelson & Zamorani (1982); SS98: 
Stevens & Strickland (1998b). 



quantify the star formation activity in galaxies. The deter- 
mination of the Sfr is fundamental to get a proper under- 
standing of the formation and evolution of the galaxies. As 
we said in the introduction, different techniques involving 
different data sets from UV to radio often yield to different 
Sfr results. Part of the problem is related with the un- 
known amount of extinction within each particular galaxy 
(|Calzetti[l200lt) , such as the amount of dust obscuration de- 
pends on the galaxy mass, galaxy type, the chemical evo- 
lutionary state, gas content or even if the galaxy is inter- 
acting or merging with another independent object. As we 
explained in the previous section, FIR and radio data pro- 
vide an extinction-free estimation of the Sfr, while FUV 
emission nicely traces the very young stellar component. 
Here, we analyzed all the available multiwavelength data 
for our sample of WR galaxies, including our reddening- 
corrected Ha estimations (see Paper I), to determine in a 
comprehensive way the Sfr within these objects. 

Table [1 compiles all FUV, U, B, Ha, H, FIR, 15 /im, 
60 /im, and 1.4 GHz luminosities for the galaxies analyzed 
in this work. The U, B, and H luminosities were com- 
puted from the reddening-corrected absolute magnitudes 
in the U, B and H bands (see Paper I) using the stan- 
dard equation logi^; — 0.4 x (M^j.fT) — M^), and consid ering 
Mu,Q=5.58 and Mr ^=5.48 (iBessell. Ga stclli fc PleJl998h 
and Mh,q=^M ([Colina. Bohlin fc Castelli.. 19961 1 

We used the values listed in Table [H to estimate the 
Sfr that each object is experiencing, following the dif- 
ferent multi-wavelength techniques explained in the previ- 
ous section. The values of the H a-based Sfr are e xtracted 
from Paper I and consider the [Kennicutn (|l998l) calibra- 
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Table 4. FUV, U, B, Ha, H, FIR, 15 fim, 60 ^m and 1.4 GHz luminosities for all galaxies analyzed in this work. 



Galaxy 


Lpuv 


Lu 


Lb 


Luc, 


Lh 


-^FIR 


Ll5fj,m 


-^60 Mm 


-^1.4 GHz 




[1039 erg s-lA- 


1] [10« Lq 


[108 Lq] 


[10*0 ergs-i] 


[10« Lq 


] [10^2 ergs-1] 


[10^2 W Hz-1 


] [1023 W Hz-1] 


[1020 W Hz-1] 


HCG 31 


49±3 


356 


202+7 


52+4 


42.2 


72+6" 


3.9+0.7" 


14.1+1.1'" 


98+14 


" AC 


24.9+1.5 


175 


92+4 


36.3+2.3 


18.0 








79+11 


" B 


11.0+0.9 


73.8 


47.4+1.8 


3.1+0.3 


10.0 








7.5+1.1 


" E 


1.51+0.12 


6.3 


3.13+0.17 


1.49+0.16 


0.586 










"F" 


1.52+0.15 


10.9 


4.0+0.2 


5.0+0.3 


0.486 










" H 


0.13+0.05 






0.045+0.013 












" G 


10.3+0.6 


90.4 


55+2 


6.4+0.4 


13.1 








11.9+1.8 


Mkn 1087 




1803 


1127+83 


70+5 


360 


228+23 


15+4 


45+5 


178+9 


" N 






26.1+0.9 


2.27+0.19 


5.01 










Haro 15 


84+5 


614 


347+13 


42+5 


124 


62+6 


11+3 


12.2+1.1 


160+9 


Mkn 1199 




479 


291+11 


49+7 


344 


116+6 


9.8+1.1 


23.8+1.2 


126+4 


" NE 




21.1 


16.6+0.6 


0.90+0.23 


21.7 










Mkn 5 


0.274+0.018 


4.21 


2.63+0.10 


0.582+0.014 


1.89 


<0.30 


<0.087 


0.036+0.007 


3.9+0.8 


IRAS 08208+2816 


162+10 


879 


511+14 


142+12 


373 


254+22 


<5.5 


49.7+4.0 


657+26 


IRAS 08339+6517 


79+4 


1159 


661+23 


120+6 


592 


198.6+1.8 


18.3+1.8 


42.6+0.3 


249+37 


" Comp. 


1.17+0.10 


38 


29.9+1.6 


2.1+0.4 


23.3 










POX 4 


7.0+0.4 


104.7 


51.1+0.5 


23.3+1.9'= 


17.9 


6.9+0.5 


<24 


1.56+0.14 


10.4+1.2 


" Comp. 


0.667+0.075 


2.42 


2.17+0.07 


0.188+0.018'= 


1.1; 










UM420 


25.4+1.4 


236 


103+4 


47+3 


57 


142: 




27.7 


74+34 


SBS 0926+606* 


41.4+3.4 


52.3 


26.6+0.8 


11.8+0.7 


9.49 


5.8+0.6 


2.8 


1.01+0.17 


10+2 


" A 


6.4+0.4 


28.1 


12.8+0.4 


9.4+0.5 


3.63 










" B 


35+3 


24.2 


13.8+0.4 


2.4+0.2 


5.86 










SBS 0948+532 




121.3 


36.7+1.0 


78+3 


15.8 








<38 


SBS 1054+365 


0.090+0.006 


0.982 


0.655+0.018 


0.450+0.017 


0.413 


0.23+0.03 


<0.042 


0.041+0.005 


0.098+0.011 


SBS 1211+540 


0.076+0.005 


0.608 


0.316+0.008 


0.141+0.005 


0.127 








<0.50 


SBS 1319+579 


0.86+0.05 


63.1 


40.2+1.4 


2.38+0.15 


28.6 


5.3: 




0.85 


<2.9 


SBS 1415+437 


0.409+0.023 


1.69 


1.00+0.03 


0.50+0.02 


0.581 








<0.052 


III Zw 107 




286 


177+5 


40.1+1.8 


77.3 


50+8 


<7.3 


10.4+1.5 


62+4 


Tol9 


9.5+0.7 


118 


79+2 


22.9+1.6 


55.5 


35.4+1.6 


2.5+0.7 


6.1+0.5 


87.3+1.6 


Tol 1457-262' 


38.6+2.6 


319 


182+5 


63+3 


109 


82+6 


<6.5 


17.2+1.1 


216+10 


" Obj 1 


29.2+1.6 


223 


121+3 


46.2+2.7 


58.6 










" Obj 2 


9.2+1.1 


97.3 


60.3+1.7 


17.2+0.8 


45.3 










" #15 


1.10+0.21 


8.39 


6.08+0.23 


0.47+0.04 


5.5: 










Arp 252'' 


47.4+3.3 


711 


435+13 


90+6 


423 


361+19 


38+5 


79+4 


1968+60 


" ESO 566-8 


38.5+2.7 


597 


350+10 


84+6 


313 










" ESO 566-7 


8.9+0.6 


114 


85+3 


6.3+0.4 


110 










NGC 5253 


2.35+0.15 


36.6 


22.9+0.2 


4.4+0.2 


14.2 


2.583+0.009 


0.479+0.004 


0.571+0.001 


1.64+0.07 



'" As IRAS data do not allow to distinguish regions within the HCG 31 group, this value considers the flux of all galaxy members. 
' We are considering the flux of two galaxies: members Fl and F2 for HCG 31 F; galaxies A and B for SBS 0926+606; obj 1 and obj 2 for 
Tol 1457-262; and ESO 566-8 (A) and ESO 566-7 (B) for Arp 252. 
As we commented in Pa per I, the Ha flux for POX 4 prov i ded b y lMendez 



EstebanI J1999I ) seems to be overestimated, hence we consider 



here the value provided bv lGil de Paz. Madore &: Pevunoval ll2003f) . that is 0.61 times smaller. We also scale our Ha flux of POX 4 Comp 
using this factor. 



tion. Recently, ICalzetti et al.' (2007) re-calibrated the rela- 
tionship between the Ha-luminosity and the Sfr; t h e Hct- 
based values of the Sfr provided bv ICalzetti et al. ( 2007 ) 
are 0.67 times the values derived using the lKennicuttI (|1998[ ) 
calibration. 

Table [5] compiles all Sfr values derived for each galaxy. 
From this table, it is evident that the agreement between 
values obtained using different methods is usually good, al- 
though sometimes we find clear discrepancies (i.e. POX 4, 
NGC 5253). We note that for systems that involve two 
or more galaxies (HCG 31, SBS 0926-^606, Tol 1457-262 
and Arp 252) we list both the global and individual Sfrs, 
because the FIR and the radio data do not have enough 
spatial resolution to distinguish the emission coming from 
different members, but FUV and Ha data do. We also con- 
sidered HCG 31 Fl and F2 as a single entity (HCG 31 F) 
because the available H I data include both TDG candi- 
dates. 

Figure [5] compares our Ha-based Sfr (corrected for 
both extinction and [N Ii] contribution as we explained in 



Appendix C of Paper I) with the Sfr estimations derived 
from FIR, 15 /zm, 60 /zm and 1.4 GHz luminosities. The 
diagram involving Li5^m seems to show a higher scatter at 
higher Sfr, but this calibration is more uncertain. As a par- 
ticular case, Arp 252 always shows a disagreement between 
the Sfr derived from Ha and other parameters, remark- 
ing with the 1.4 GHz luminosity. The main object within 
Arp 252 is the bright galaxy ESO 566-8. This behavior, to- 
gether the fact that the FZR-radio-continuum relation is not 
satisfied in this system (see below) strongly suggest that 
ESO 566-8 has some activity different to its starbursting 
nature (an AGN or a radio-galaxy), something we already 
commented when we analyzed this system (see Sect. 3.19.2 
of Paper I). The rest of the objects agree rather well when 
comparing values obtained from different calibrations. As 
previous authors pointed out (i.e. Dopita et al. 2002; James 
et al. 2005), the correction of the Ha fiuxes for both extinc- 
tion and [N ii] emission is vital to get a reliable estimation 
of the SFR using Ha-images. 
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Tabl e 5. Sfr values (in units of Mp yr ^ ) derived for ea ch galax y using different lu minositi es and calibrat i ons. T he refere nces are: 
S07 =ISaIInreraII ||2007|); G84 =lGa'llagher et al.1 (Il98^ ): K98 = lKennicuttl lfl99g ): ROl = iRoussel et"aLl (|200lh ; C92 = lCondonl 
l|1992r ) ; C02 = ICondon et alJ (|2002l ) using the expression for M >5 Mq. The last column indicates the assumed value of the Sfr 
we estimate to each galaxy considering all available data. 



Galaxy 



FUV 



U 



B 



Ha 



FIR 



1.4 GHz 



Assumed 





S07 


Eq.[T5] 


G84 


Eq.IIl 


K98" 


K98 


ROl 


C92 


C02 


Sfr'' 


HOG 31 


4.0+0.3 


2.8 


0.59 


2.5 


4.1+0.4 


3.3+0.3 


1.2+0.2 


2.8+0.2 


2.5+0.2 


3.1 


" AC (NGC 1741) 


2.02+0.12 


1.5 


0.27 


1.2 


2.88+0.18 








1.98+0.03 


2.0 


" B 


0.89+0.07 


0.68 


0.14 


0.67 


0.24+0.02 








0.19+0.03 


0.62 


" E 


0.122+0.009 


0.073 


0.009 


0.057 


0.118+0.013 










0.10 


" F 


0.123+0.013 


0.12 


0.012 


0.071 


0.40+0.03 










0.20 


" G (Mkn 1090) 


0.83+0.05 


0.81 


0.16 


0.77 


0.51+0.03 








0.30+0.04 


0.49 


" H 


0.011+0.004 








0.004+0.001 










0.008 


Mkn 1087 




12.1 


3.27 


11.8 


5.6+0.4 


10.3+1.1 


4.6+1.3 


8.8+1.0 


4.5+0.2 


6.3 


" N 






0.076 


0.39 


0.180+0.015 










0.12 


Haro 15 


6.8+0.4 


4.6 


1.00 


4.1 


3.3+0.4 


2.8+0.3 


3.2+0.9 


2.4+0.2 


4.0+0.2 


3.6 


Mkn 1199 




3.7 


0.84 


3.5 


3.9+0.6 


5.2+0.3 


3.0+0.3 


4.7+0.2 


3.16+0.10 


3.7 


" NE 




0.22 


0.048 


0.26 


0.07+0.02 










0.05 


Mkn 5 


0.0222+0.0014 


0.050 


0.008 


0.049 


0.046+0.011 


<0.014 


<0.03 


0.067+0.008 


<0.10 


0.040 


IRAS 08208+2816 


13.1+0.8 


6.4 


1.48 


5.8 


11.3+0.9 


11.4+1.0 


<1.68 


9.7+0.8 


16.4+0.7 


11.6 


IRAS 08339+6517 


6.4+0.4 


8.2 


1.92 


7.3 


9.5+0.5 


8.93+0.08 


5.6 + 0.6 


8.36 ± 0.06 


6.2+0.9 


7.0 


" Comp 


0.095+0.008 


0.37 


0.087 


0.44 


0.17+0.02 










0.10 


POX 4 


0.57+0.03 


0.93 


0.15 


0.72 


1.85+0.06 


0.31+0.02 


<7.48 


0.31+0.03 


0.26+0.03 


0.54 


" Comp 


0.054+0.006 


0.030 


0.0063 


0.041 


0.012+0.004 










0.031 


UM420 


2.01+0.11 


1.9 


0.30 


1.4 


3.7+0.2 


6.4<= 




5.4^= 


1.9+0.8 


2.1 


BBS 0926+606 


3.4+0.3 


0.50 


0.08 


0.40 


0.94+0.06 


0.26+0.03 


<0.86 


0.20+0.03 


0.25+0.06 


0.95 


" A 


0.53+0.03 


0.28 


0.04 


0.20 


0.75+0.04 










0.52 


" B"* 


2.82+0.25 


0.25 


0.04 


0.22 


0.19+0.02 










1.4? 


SBS 0948+532 




1.1 


0.11 


0.53 


6.2+0.2 








<0.95 


4.2 


SBS 1054+365 


0.0073+0.0005 


0.013 


0.0014 


0.019 


0.036+0.001 


0.016 


<0.01 


0.015 


0.025+0.004 


0.018 


SBS 1211+540 


0.0062+0.0004 


0.009 


0.0007 


0.009 


0.011+0.001 








<0.01 


0.007 


SBS 1319+579 


0.069+0.004 


0.59 


0.12 


0.58 


0.189+0.012 


0.24 




0.17 


<0.07 


0.15 


SBS 1415+437 


0.0331+0.0018 


0.022 


0.0029 


0.020 


0.039+0.002 








<0.01 


0.030 


III Zw 107 




2.3 


0.51 


2.2 


3.19+0.15 


2.3+0.4 


<2.24 


2.0+0.3 


1.52+0.09 


2.0 


ToI9 


0.77+0.06 


1.0 


0.23 


1.1 


1.82+0.13 


1.59+0.07 


0.8+0.2 


1.19+0.08 


2.18+0.04 


1.3 


Tol 1457-262 


3.2+0.2 


2.6 


0.53 


2.3 


5.0+0.3 


3.7+0.3 


<1.99 


3.4+0.2 


5.4+0.3 


3.8 


" Obj 1 


2.37+0.13 


1.8 


0.35 


1.6 


3.7+0.2 










2.4 


" Obj 2 


0.74+0.09 


0.87 


0.17 


0.83 


1.37+0.06 










0.83 


" #15 


0.089+0.017 


0.09 


0.018 


0.10 


0.038+0.003 










0.06 


Arp 252 


3.84+0.26 


5.3 


1.26 


5.0 


7.2+0.5 


16.2+0.9 


11.6+1.4 


15.4+0.8 


49.2+1.5 


10 


" ESO 566-8 


3.12+0.22 


4.5 


1.01 


4.1 


6.7+0.5 










3.8 


" ESO 566-7 


0.72+0.04 


1.0 


0.25 


1.1 


0.50+0.04 










0.53 


NGC 5253 


0.190+0.010 


0.36 


0.07 


0.35 


0.348+0.017 


0.12 


0.15 


0.11 


0.041+0.002 


0.14 



" Considering the new correlation between the Sfr and the Ha-luminosity provided bv lCalzetti et al] ||2007| ). the Ha-based Sfrs arc 0.67 

times the values shown here. 

' We consider the Ha-based Sfr. values provided by the lCalzetti et aU ||2007I ') calibration to estimate the average Sfr (see text). 

'^ The FIR and 60 /xm luminosities in UM 420 are overestimated because of the contribution of the foreground galaxy UGC 01809 (see 

Fig. 13 and Sect. 3.9.1 in Paper I). 

■^ SBS 0926+606 B shows extended FUV emission, as it is seen in the right panel of Figure^ The values of the Sfr provided by its FUV 

and Ha emission are very different. 



Although the agreement between the Ha-based Sfr and 
the Sfrs derived using FIR and radio luminosities is good, 
we observe that the values provided using the Ha luminos- 
ity are slightly higher than those estimated using the other 
calibrations. The difference seems to be higher at lower Ha- 
luminosities. A linear fit to the data (green continuous lines 
in Figure ^ confirms this trend. The zero-points of the fits 
(0.59, 0.52, 0.58, 0.68 for the Ra-FIR, Ha-1.4 GHz, Ha- 
60 /um, and Hq;-15 fim relations, respectively) indicate that, 
for Sfr=1 Mq yr^^, the value of the Sfr provided by Ha- 
luminosity is ~0.6 times the Sfr values estimated using the 
other relations. Bell (2003) concluded that both radio and 
FIR luminosities underestimate the Sfr for low-luminosity 
galaxies because the non-thermal emission seems to be sup- 
pressed by a factor of 2-3 in dwarf objects. However, the 
difference is not significative if we use the ICalzetti et al.l 



()2007D calibration instead of the iKennicuttJ ()1998f ) calibra- 
tion to derive the Ha-based Sfr. 

The comparison of the FUV-hased with the Ha-based 
Sfr (Figure [3]) also shows a good agreement: except for 
some few objects (remarkably SBS 0926+606 eQ) both 
relations provide similar values. We also observe that 
the FUV-hased Sfrs seem to be slightly lower than 
the Ha-based Sfrs. A linear fit to the data (shown in 
Figure [3] with a continuous green line and with a cor- 
relation coefficient of r=0.927) indicates that the FUV- 
hased Sfr is, on average, '^0.71 times the Ha-based Sfr. 
This value is similar to the factors found before when 
comparing the Ha-based Sfr with the FIR- and radio- 



^ As we commented before, the FUV emission observed in 
SBS 0926+606 B is much more extended that the Ha emission, 
and hence the derived Sfr is more than one order of magnitude 
higher using the FUV than the Ha emission. 
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Fig. 2 . Comparison between the Ha-based Sfr -corrected for both extinction and [Nil] contribution and assuming iKennicutt] 
l|l998t ). X-axis- with the Sfrs derived using the FIR, 15 /^m, 60 /im and 1.4 GHz luminosities. The dotted red lines indicate the 
position with equal Sfr, the continuous green lines show a lineal fit to the data. 




based Sfrs. Interestingly, all the se numbers are c oinci- 
dent with the ratio between the iKennicuttI (I1998D and 



SFR (Ha) 

Fig. 3. Comparison between the Ha-based Sfr -co rrected for 
both e xtinction and [Nil] contribution and assuming iKennicutti 
i|1998t) . x-axis- with the Sfrs derived using the FUV luminosi- 
ties. The dotted red line indicates the position with equal Sfr, 
while the continuous green line shows a fit to the data. 



the Calz etti et ahl ()2007[ ) calibrations to the Sfr using the 
Hct flux, S'Fi?co7(Ha)/S'Fi?K98(Ha)=0.67. We therefore 
conclude that the n ew Ha-based calibration provided by 
ICalz etti et al.' ('2007') should be preferred over the widely- 
used Kcnnicutt (1998) calibration when computing the Sfr 
from Hta luminosities. The Sfr estimated for each object 
considering all available multiwavelength data, and listed in 
last column of Table O has been computed considering the 
ICalzetti et al.l (|2007[ ) value. Finally, we must say that there 
are increasing evidences that the Ha luminosity underes- 
timates the Sfr relative to the FUV luminosity in dwarf 
galaxies with Sfr<0.01 Moyr"^ (i.e., Lee et al. 2009) and 
hence the Ft/F-based Sfr should be preferred over the 
Ha-based Sfr in those systems. 

3.1. Ls-Sfr a/io' L^z-Sfr relations for starburst galaxies 

Just for comparison, we also estimated the Sfr from the B- 
luminosity using the calibration provided by Gallagher et 
al. (1984). SFRb represents the star formation activity oc- 
curred in the last few hundreds of Myr, while the rest of the 
calibrations are tracing the massive stars and the nebular 
emission of the gas that only last for some few tens of Myr. 
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Fig. 4. Assumed Sfr vs. B-luminosity (left panel) and [/-luminosity (right panel) for the analyzed galaxies. Luminosities are 
plotted in solar units. The best fit (in logarithm scale) to our data are plotted with a continuous red line. The previous calibration 
given by Gallagher et al. (1984) between the Sfr and the _B-luminosity is shown by a discontinuous green line. 
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Fig. 5. Sfr vs 12-|-log(0/H) {left) and Sfr/-Ls vs 12-|-log(0/H) (right) for our sample of WR galaxies. The red-dotted line 
indicates a fit to our data. Green diamonds in left panel plot the average value obtained in the low, intermediate and high- 
metallicity regimes. 



For our galaxy sample, SFRb is always lower than the Sfr 
derived from the other calibrations, as we should expect be- 
cause of the starbursting nature of the analyzed galaxies. 
The value of the Sfr^ in Mkn 1087 using Gallagher et al. 
(1984) equation is only half of that estimated from other 
calibrators, remarking its Luminous Blue Co mpact Galaxy 
(LCBG) nature (Lopez-Sanchez et al.|[2004bl) . 



We used our data to establish a new relation between 
the Sfr and the -B-luminosity, that should be applied only 
in starburst galaxies and just as a first estimation of the 
actual Sfr. The left panel of Figure S] shows the rela- 
tion between Lb (in solar units) and the assumed Sfr for 
all our galaxies. Despite some clear discrepancies between 
some galaxies that show very different Sfr for a similar B- 
luminosity (for example, just compare members G and F 
of HCG 31), we observe a good agreement, having galaxies 
with higher -B-luminosities higher star-formation activity. 
The discrepancies are consequence of the different star- 
formation histories of the galaxies (relative contribution 
and age of the underlying stellar population, metallicity. 



age of the most recent star- formation event). A linear fit to 
our data provides the relation 



S'-F-R_B starbursts — 1.148 X 10 L 



-9r0.906 



(14) 



being Lb expressed in units of Lq. The correlation coef- 
ficient of this fit is r=0.932. The third column in Table [5] 
compiles the SFRb computed for each galaxy using this 
equation. The relation obtained by Gallagher et al. (1984) 
gives values one order of magnitude lower than those we 
obtain with our new calibration. 

We also computed a relation between the Sfr and the 
L'^-luminosity for our sample galaxies. The right panel of 
Figure m shows such relation. A linear fit to the data yields 



SFR 



t/,starbursts 



7.59 X 10-^°L' 



0.907 
U ' 



(15) 



being Ljj expressed in units of Lq. This fit has somewhat 
smaller scatter than our derived Sfr-Lb relation, resulting 
in a correlation coefficient r=0.965. However, the slope for 
both calibrations (0.894±0.070 and 0.900±0.054 for Sfr- 
Lb and Sfr-Lij, respectively) are similar. We remark that 
Equations [T3] and [1^] can not be applied to galaxies with 
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no strong star-formation activity, because the Sfr value 
derived from them will be overestimated (as it is happening 
in Mkn 1199 NE or IRAS 08339-1-6517 comp). 



3.2. Comparison of Sfr and metallicity 

Left panel of Figure [5] compares the assumed Sfr with the 
oxygen abundance computed for each galaxy. We estimated 
the average Sfr values in the low (12-|-log(0/H)<7.8), 
intermediate (7.8<12-Mog(0/H)<8.3) and high 
(12-t-log(0/H)>8.3) metallicity regimes, the we plot 
in this figure using green diamonds. As we see, the 
dispersion in the intermediate-metallicity range is quite 
high, but that is just a consequence of the star-formation 
history of each particular galaxy (see Paper IV), as in this 
metallicity regime lies both very dwarf objects (i.e., Mkn 5, 
SBS 1054-f 365) and large and bright star-forming galaxies 
(i.e., Tol 1456-262, III Zw 107) which share a relatively 
similar chemical history. Besides the large dispersion in the 
intermediate-metallicity regime, it is clear that galaxies 
with higher metallicity have higher global star formation 
rates. That is a consequence of the building of the galaxies, 
because more massive objects are more metal-rich than 
less massive galaxies (see below) and, hence, when the 
starburst is initiated, galaxies with higher mass (and with 
higher metallicities) will create stars at a higher rate 
than those found in smaller objects. The comparison of 
the Sfr per B-luminosity, Sfr/Lb with the metallicity 
(Figure [5l right) also shows a tremendous dispersion for 
12-|-log(0/H) between 8.0 and 8.2. However, we observe 
that Sfr/Lb decreases with increasing oxygen abundance 
(the red-dotted line shows a fit to the data), indicating 
that galaxies with lower metallicity (and, therefore, less 
massive objects) have stronger star-forming bursts than 
those found in higher metallicity (more massive) objects. 
SBS 0948-1-532 has the highest Sfr/Lb in our sample, 
indicating the strength of the starburst, as we saw when 
analyzed its photometric properties (see Sect. 3.11 in 
Paper I). On the other hand, SBS 1319+579 has a very 
low Sfr/Lb in comparison with BCDGs of similar char- 
acteristics, indicating the peculiarity of this galaxy. We 
will see below that other properties of SBS 1319+579 
show additional discrepancies with the average behavior 
in BCDGs, suggesting that the star-formation activity has 
been somewhat suppressed in this object. For example, the 
gas depletion timescale is extremely long for a starburst 
galaxy, (r ~12.7 Gyr, see Table IH]). 



3.3. A Lx-Sfr relation for starburst galaxies 

Although several relations between the X-ray luminos- 
ity and the star formation rate have been proposed (i.e., 
Ranalli et al. 2003; Lou & Bian 2005) they seem not 
to be appropriate for young starbursting systems. For 
example, as we ex plained in the analysis of the LCBG 
IRAS 08339+651 7 (iLopez-Sanchez et all l2006l) . the rela- 
tion provided by iRanalli et al.l (|2003l) gives a very high 
Sfr value (61.8 M0yr~^) in comparison with the esti- 
mations obta ined using other frequencies (6-8 M0yr~^). 
IStevens fc Str ickland ( 1998a ) showed that the X-ray lumi- 
nosities in WR galaxies are substantially higher that those 
found in non-WR galaxies with similar i?-luminosity. That 
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Fig. 6. X-ray luminosity in the 0.2-2.0 keV range vs. FIR lu- 
minosity for the sample of WR galaxies analyzed by Stevens & 
Strickland (1998a,b). The red continuous line is the best fit to 
the data, excluding the values for NGC 5253 and NGC 5408. 
The green discontinuous line is th e relation ob tained using 
the Sfr-Lx calibration provided bv iLou fc BiarJ ((2005), while 
the yellow dotte d -dash ed line is the relation obtained from 
the IRanalli et al.l (|2003r i calibration. The three additional WR 
galaxies of our sample for which X-ray data are available 
(HCG 31, IRAS 08339+6517, and Tol 9) are indicated with dark 
yellow squares. 



is a consequence of the higher rate of superbubbles and 
supernova explosions in WR galaxies. 

We have used the sample of WR galaxies analyzed by 
Stevens & Strickland (1998a, b) to get a tentative calibra- 
tion between Sfr and Lx for this kind of objects. These 
authors obtained X-ray data in the 0.2-2.0 keV range us- 
ing the satellite ROSAT. We have checked which of these 
galaxies also possess FIR data from the IRAS satellite, 
and established a relation between Lpm and Lx, as it is 
shown in Figure |6l Only 18 galaxies have available data for 
both luminosities. NGC 5253 was included in the Stevens & 
Strickland (1998a, b) analysis, but they indicated that the 
X-ray emission in this object is very peculiar. The X-ray 
emission measured in NGC 5408 may be unrelated with 
the galaxy. Hence, neglecting the contribution of these two 
galaxies, the linear fit to the data gives 



L 



FIR 



= 24x L 



X 



(16) 



being the correlation coefficient r=0.9 29. Considering the 
calibration given by Kcnnicuti ()1998f) between LpiR and 
Sfr (EquationH]), we find the following calibration between 
soft X-rays (0.2-2.0 keV) and Sfr: 



SFRx = 1.08 X 10 



-42 7-1.04 
^X ■ 



(17) 



As it can be seen in Figure [6l our new Sfr-Lx calibration 
agrees better with the ob servational data tha t the relations 
explained before. Indeed, iLou fc BianI (|2005[ ) relation gives 
values around one order of magnitude lower than those ex- 
pected from the FIR luminosity, but IRanalli et al.l ()2003[) 
calibration provides values almost one order of magnitude 
higher than the actual ones. Using our new relation and the 
available X-ray data for our WR galaxies (see Table ^ , we 
derive a SFRx of 3.5, 13.8, and 0.86 Mq yr"! for HCG 31, 
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Table 6. Additional FIR and radio properties. 



Galaxy 


<?" 


logBt' 


[Gyr] 


HCG 31 AC 


2.39 


1.19 


1.92 


HCG 31 B 




1.25 


5.3 


HCG 31 F 






2.9 


HCG 31 G 




1.11 


2.3 


Mkn 1087 


2.53 


1.30 


3.3 


Haro 15 


2.01 


1.48 


1.9 


Mkn 1199 


2.39 


1.31 


0.40 


Mkn 5 






2.1 


IRAS 08208+2816 


2.01 


1.57 




IRAS 08339+6517 


2.33 


1.21 


0.89 


POX 4 


2.25 


0.93 


4.2 


UM420 


2.71 


1.08 




SBS 0926+606 


2.18 


0.79 


2.0 (A), 5.6 (B) 


SBS 0948+532 








SBS 1054+365 


2.79 


-0.10 


4.0 


SBS 1211+540 






3.5 


SBS 1319+579 


2.15 




12.7 


SBS 1415+437 






3.5 


HI Zw 107 


2.34 


1.06 


3.5 


TolQ 


2.03 


1.16 


1.95 


Tol 1457-262 


2.00 


1.43 


1.51 


Arp 252 


1.69 


2.28 




NGC 5253 


2.62 


0.33 


0.97 



'^ The logarithmic ratio of FIR to radio flux density parameter, q, i 

defined in Equation 1201 

^ R'ls the non-thermal to thermal ratio, derived from the 1.4 GHz 

and Ha fluxes. 

■^ The gas depletion timescale is defined as Tg = 1.32A/hi/Sfr 

(Skillman et al. 2003) 



IRAS 08339+6517, and Tol 9, respectively. These values 
agree well with the actual Sfr estimated for each object 
(see Table [5]) . 



4. FIR/radio correlation 

We have used the luminosity data shown in Table|3]to check 
if our WR galaxies follow the FIR/r&Aio correlation. As it 
was shown by Condon et al. (1992), the FZR/radio correla- 
tion is much tighter for starbursts than for active galaxies. 
Figure [7] (left) plots the 1.4 GHz luminosity vs. the 60 /im 
luminosity for our sample galaxies and the relation betw een 
both quantities found bv lYun. Reddv fc CondonI (|200l[ ). 



ii.40GHz [W Hz-1] - lO^^Leo ^™ [L©], (18) 

while Figure [7] (right) shows L-\An \^v. vs. the ioia l FIR hx- 
minosity and the relation given by ICondon et alj (|1991f ) , 



log Li.49GHz[W Hz-1] = 1.1 log Lfir [Lq] + 10.45. (19) 

Bell (2003) remarked that the radio- f/i? correlation is lin- 
ear not because both radio and FIR emission track Sfr, 
but rather because they fail to track Sfr in independent, 
but coincidentally quite similar, ways. Further analysis (i.e., 
Hunt et al. 2005) also found that this relation is not hold for 
some low-metallicity or young starbursts galaxies. However, 
as it is seen in Figure [71 all analyzed objects except Mkn 5 
(which has a very uncertain value for FIR) and Arp 252 
(ESO 566-8 hosts some kind of nuclear activity) follow well 
both relations. This indicates that the galaxies are star- 
bursting systems and are not active galaxies (Seyfert or 
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Fig. 8. Comparison of the i3- luminosity (in solar units) and the 
logarithmic non-thermal to thermal ratio, log -R, for our sample 
galaxies. A linear fit is shown with a continuous red line. 



AGNs). We already reached this conclusion when we an- 
alyzed the diagnostic diagrams involving several emission- 
line ratios (see Paper HI). Figure [7] includes a linear fit (in 
logarithmic scale) to our data (neglecting Mkn 5, for which 
th e FIR values have hig h uncertainties) . The relation given 
bv lCondon et al.1 (|1991[ ) seems to be slightly displaced with 
respect our observationa l data, although we also see so me 
small discrepancies in the lYun. Reddv fc GondorJ (J200lf l re- 
lation for the faintest objects. 

The non-AGN nature of our sample of WR galaxies is 
also supported by the analysis of the q parameter an the 
FIR spectral index. The q parameter is defined as the log- 
arithmic ratio of FIR to radio flux density. 



lot 



FfibXW m-2) / 3.75 x lO^^ Hz 
^1.4 GHz (W m-2 Hz-i) ' 



(20) 



and it is very robust for most galaxy populations; < g > = 
2.34 ± 0.19 (Condon et al. 1991; Yun et al. 2001). Galaxies 
with q < 1.8 are more than 3 times more radio loud than the 
mean for the star-forming galaxies, so they can be classified 
as AGN-powered. As it seen in Tabled all galaxies except 
Arp 252 have a q value similar to that expected for starburst 
galaxies. 

Table ini also compiles the non-thermal to thermal ratio, 
R, of the galaxies with available 1.4 GHz radio-continuum 
data. The thermal flux at 1.4 GHz was computed apply- 
ing Equation [5] The majority of the galaxies show the 
typical value for st ar-forming galaxies, \ogR = 1.3 ± 0.4 
iDopita et all (|2002D . The low value in R found in POX 4 
and NGC 5253 may be because the Ha flux has been 
overestimated, although the situation of NGC 5253 is far 
from clear (Lopez-Sanchez et al. 2010a). The value obtained 
for SBS 1054+365 is not reliable, we consider that or the 
1.4 GHz flux was underestimated (very probably) or the 
Ha flux was overestimated. However, the high value found 
in Arp 252 (the emission comes mainly from ESO 566-8), 
log-R = 2.28, is real and indicates that the thermal flux at 
1.4 GHz is less than 0.5%. As reported by several authors 
(i.e., Klein, Wielebinski & Thuan 1984; Klein, Weiland & 
Brinks 1991; Bell 2003), dwarf galaxies seem to have a lower 
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Fi g. 7. 1.4 GHz radio-con t inuum luminosity vs. the 60/im lumino sity (left) and the FIR luminosity (right). The relations derived 
bv lYun. Reddv fc Condon! (|2001f ) (Equation 1181 left diagram) and lCondon et al] (| 19911 ) (Equation 1191 right diagram) are plotted 
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non-thermal-to-thermal emission ratio than normal spiral 
galaxies. The values obtained for the R parameter in our 
galaxy sample tend to be lower at lower B-luminosities, 
as it is shown in Figure [S] The difference between dwarf 
and larger galaxies is often interpreted as higher efficiency 
of cosmic-ray confinement in more massive galaxies (e.g., 
Klein et al. 1984; Price & Duric 1992; Niklas, Klein & 
Wielebinski 1997; Beh 2003). 



5. Analysis of the masses 

For this work, we have estimated the ionized gas mass - 
-A^Hii) using the Ha images presented in Paper I-; neutral 
gas mass --Mhi, using Hi data at 21 cm compiled from the 
literature-; mass of the ionizing star cluster -M^, using Ha 
and VF(H/3), see Paper I-; warm dust mass --Mdust, using 
the FIR fiuxes-; Keplerian mass -./V/kop, via the kinematics 
of the ionized gas-; and dynamical mass -.Mdyn, using the 
H I kinematics-. All these data are compiled in Table[71 The 
estimation of -Mxep and -Mdyn for each galaxy was explained 
in Paper H. We just remember that, as the extension of 
the neutral gas is usually larger than the stellar compo- 
nent, our estimations of -Mdyn are very probably underesti- 
mated. Furthermore, no-rotational movements would yield 
in a overestimation of the total mass. Only interferometric 
H I analysis can definitely provide a more precise determi- 
nation of the dynamical mass of each system. However, we 
may use our Mdyn values as a rough estimation of the total 
mass of the systems. Their comparison with Afxcp, -^^hi, 
-Mhii, -Mdust and their associated mass-to-light ratios will 
give clues about the galaxy type, dynamics and the fate of 
the neutral gas. 

We first compare all mass determinations with the opti- 
cal luminosity of the galaxies. Figure IHl shows the relations 
between Mhii, Mhi, -Mdust and Mi, with the absolute B- 
magnitude. As we should expect, and besides some scatter, 
all mass determinations clearly increase with increasing op- 



tical luminosity. We have performed a linear fit to the data; 
the results are: 

log Mhii = (-0.21 ± 0.56) - (0.342 ± 0.030)Ms, (21) 

log Mhi = (3.69 ±0.52)- (0.292 ±0.G28)Ms, (22) 

log Mdust = (-1.15 ± 0.92) - (0.364 ± 0.047)Ms, (23) 

log M* = (-0.47 ± 0.54) - (0.397 ± 0.029)Mb, (24) 

with correlation coefficients of 0.899, 0.922, 0.912 and 
0.928, respectively. Some deviations to the fits are found 
in Mkn 1199 (that possesses a relatively low Mhi), 
SBS 0948-^532 (with a very high Afnii), POX 4 (it seems 
to be Mdust deficient, while its ionized gas mass may be 
overestimated), UM 420 (its high Mdust is very proba- 
bly consequence of the contamination of the FIR emis- 
sion by the foreground galaxy UGC 01809, see Sect. 3.9 
in Paper I and Sect. 3.9 in Paper II), SBS 1319+579 and 
IRAS 08339+6517 Comp (that have a very low A/* for their 
absolute i?- magnitude) and NGC 5253 (that is both Afni 
and Afdust deficient). 

Figure [TU] plots the dynamical mass (that represents the 
total mass of the galaxy) versus the absolute magnitude in 
several broad-band filters (.B, F, R and J). We find a clear 
correlation between these quantities, a linear fit to the data 
yields: 



log Mrfj,„ = (4.37 ±0.58) - 
log Md„„ = (4.28 ±0.55) - 
log M<i„„ = (3.86 ±0.57) - 
logAf<ij,„ = (4.60 ±0.65) 



(0.304 ±0.032)Mb, (25) 

(0.306 ±0.030)My, (26) 

(0.324 ± 0.030)Mi^, (27) 

- (0.281 ± 0.034)Mj, (28) 



with correlation coefficients r of 0.922, 0.931, 0.940 and 
0.907, respectively. Notice that slopes in all fits are 
quite similar. The most important deviations to these fits 
are found in clearly interacting systems (Mkn 1199 and 
HCG 31 AC) but also in Mkn 5 and SBS 1054+365. 

We have compared the Keplerian mass (derived from 
the kinematics of the ionized gas) with the dynamical 
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Table 7. Keplerian mass (Mkop), dynamical mass (Mdyn), neutral gas mass (Afni), ionized gas mass (Mhii), warm dust mass 
(Mdust), mass of the ionizing star cluster (Af*), total stellar mass {Mstars), and baryonic mass (Mbar) of the galaxies analyzed in 
this work. 



Galaxy 


AfKcp 


A/dyn 


Mhi 


Mhii 


A/dust" 


M-, 


Mstars 


A4ar 




[10** Mq] 


[10« Mq] 


[10« Mq] 


[10'' Mq] 


[106 Mq] 


[10« Mq] 


[10« Mq] 


[10» Mq] 


HCG31 






74.6 


7.8+0.5 


2.56 


37.8+2.8 


33.8 


141 


" AC 


340 


850 


36.4 


5.4+0.3 




18.2+1.2 


UA 


62.5 


" B 


26 


54 


19.4 


0.45+0.04 




8.2+0.7 


8.0 


33.6 


" E 




13'' 




0.22+0.02 




2.6+0.3 


0.47 




" F 


3.0 


15 


6.13 


0.75+0.04 




0.40+0.03 


0.39 


8.5 


" G 


21 


68 


19.4 


0.96+0.05 




7.0+0.4 


10.5 


36.1 


Mkn 1087 


560 


1800 


156 


10.4+9.7 


7.8 


129+9 


288 


494 


" N 


2.2 


45'' 




0.337+0.028 




3.34+0.28 


4.0 




Haro 15 


121 


365 


55±18 


6.2+0.7 


2.1 


77+9 


99.2 


172 


Mkn 1199 


82 


1650 


12.2 


7.3+1.1 


3.1 


140+21 


275 


291 


" NE 


2.9 


70'' 




0.13+0.03 




1.7+0.4 


17.4 




Mkn 5 


21 


36 


0.72+0.09 


0.081+0.010 


0.099 


0.300+0.012 


1.5 


2.5 


IRAS 08208+2816 


39 


600** 




21.1+1.7 


8.84 


166+14 


298 




IRAS 08339+6517 


100 


370 


53+6 


14.5+0.7 


3.91 


226+10 


476 


546 


" Comp 


80 


100 


7.0+0.9 


0.26+0.03 




4.1+0.5 


18.6 


27.9 


POX 4 


5.0 


76 


21 


5.70+0.18 


0.093 


9.8+0.3 


14.3 


29.1 


" Comp 








0.05 




0.79+0.07 


0.88 




UM 420 


21 


200** 




6.9+0.4 


5.0^= 


13.8+0.9 


45.6 




SBS 0926+606 






17.7+7.2 


1.75+0.10 


0.37 


8.4+0.6 


7.59 


30.9 


" A 




23 


9.6+3.6 


1.40+0.07 




3.59+0.17 


2.9 


15.5 


" B 




45 


8.1+3.6 


0.35+0.03 




4.8+0.4 


4.7 


15.4 


SBS 0948+532 


21 


90'' 




11.6+0.4 




18.8+0.7 


12.6'' 




SBS 1054+365 


0.78 


15 


0.61+0.06 


0.067+0.003 


0.012 


0.23+0.01 


0.33 


1.13 


SBS 1211+540 


1.13 


1.14 


0.24+0.04 


0.021+0.001 




0.05 


0.102^* 


0.420 


SBS 1319+579 


86 


140 


16.4 


0.35+0.02 


0.30 


0.78+0.05 


22.9 


44.5 


SBS 1415+437 


2.5 


4.9 


0.96+0.07 


0.074+0.003 




0.14+0.01 


0.46 


1.74 


III Zw 107 


8.2 


180 


67+12 


6.0+0.3 


1.25 


74+3 


61.8 


150 


Tol9 


12 


580*^ 


22+2 


3.4+0.2 


2.41 


27.4+1.9 


44.4 


106 


Tol 1457-262 




950 


47 


9.4+0.5 


1.83 


41+2 


87.2 


149 


" Obj 1 


62 


290'' 




6.9+0.4 




23.4+1.5 


46.9 




" Obj 2 




150** 




2.56+0.12 




14.4+0.7 


36.2 




Arp 252 








13.4+0.9 


6.32 


82.3+4.9 


338 




" ESO 566-8 


73 


440'' 




12.5+0.8 




45+3 


250 




" ESO 566-7 


4 


150* 




0.93+0.07 




16.8+1.2 


88 




NGC 5253 




83'= 


1.63+0.10 


0.65+0.03 


0.042 


2.22+0.11 


11.4 


13.6 



" This value is for the entire system: all galaxies in the HCG 31 group, members A and B in SBS 0926+606, all galaxies in Tol 1457-262, and 

ESO 566-8 and ESO 566-7 in Arp 252. We neglect the contribution of the FIR emission in dwarf objects associated to larger galaxies 

(companion objects surrounding Mkn 1087, Mkn 1199, IRAS 08339+6519, and POX 4). 

'' Tentative value of Afdyn computed using Equations I25I428I and 1291 

'^ The warm dust mass is very probably overestimated because of the FIR contribution of the foreground galaxy UGC 01809. 

"^ Mstars computed assuming V — J ^ 0.8 and J — H ^ 0.3. 

"^ The detailed analysis of the H I gas and its kinematics will be presented elsewhere (LS+10). 



mass (estimated from the kinematics of the neutral gas). 
Figure [TT] plots both sets of values. As we expected, M^cp 
is lower than Mdyn for almost all cases {MKep^M^yn is 
shown by a dotted green line in Figure [TT]). Although the 
dispersion is high -and we remember that A^Kep and/or 
Mdyn may be overestimate because of interaction features- 
we have performed a tentative fit to the data, that yield 



Mdyn = 0.584 X M 



1.114 
Kep • 



(29) 



with a correlation coefficient r — 0.827. This relation is 
included in Figure [TT] as a red dashed line. As we explained 
in Sect. 3.13 of Paper II, Mxcp in SBS 1211-1-540 has been 
probably overestimated, so we did not consider this point in 
the analysis. The relation indicates that Afxcp is between 

12% (for Afdyn = 10" Mq) and 24% (for A/dyn = 10** Mq) 

the total dynamical mass. Hence, as we should expect, more 
massive galaxies have a higher Mdyn/MKop ratio, indicating 
that the kinematics of the ionized gas is not appropriate to 
derive the total dynamical mass in those objects. 



Using Eg nations [25ti28l and 1291 we have computed a ten- 
tative value for the dynamical mass in the galaxies with lack 
of H I data. We have included the results in Table [71 and 
plotted these points in Figure [TT] (dark yellow diamonds) . 
As we can see, they match well with the positions of the 
galaxies for with Mdyn was derived from H i data, but we 
will not consider these points in the subsequent analysis. 

We prefer to use our NIR data to derive a proper value 
of the stellar mass o f all the galaxies Following the de- 
scription provided bv iKirbv et al.l (|2008[ ). we may assume 
a il-band mass-to-light ratio of Mgtars/LH— 0-8 to com- 
pute the stellar mass, Mstars, from the _ff-luminosity (that 
is compiled for all objects in Table [p. This assumption is 
well s upported by both observatio ns (|Bellll2003tlKirbv et al.l 
200^ and theory (jde Jond 11996] ), and considers a 12 Gyr 
old solar metallicity stellar population with a constant star- 
formation rate and Salpeter initial mass function. Hence, 
the -ff-band mass-to-light ratio may be somewhat overesti- 
mated for our young galaxies. Combining the -ff-band de- 
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Fig. 9. Ionized gas mass (Mhu), neutral gas mass (Mhi), mass of the ionizing star cluster (M*), and warm dust mass (Mdust) 
vs. the absolute B magnitude for the analyzed galaxies. Linear fits to the data are shown with a dashed red line. 



rived stellar mass and the H i mass (we neglect the ionized 
gas, molecular gas, and dust contributions), the total bary- 
onic mass, Mbar, can be computed via 



Mbar — Ms 



1.32A/HI 



(30) 



where the factor 1.32 corrects the H i mass for the presence 
of hehum. The derived values for both Aibar and M^tars are 
compiled in last columns in Table |7l For SBS 0948-1-532 and 
SBS 1211-1-540, that lack of NIR colors, we have assumed 
that y - J '- 0.8 and J - H ^ Q.ito derive the 7J-band 
luminosity. 

As we should expect, the comparison between the 
dynamical and the baryonic masses (Figure I12p indi- 
cates that Mdyn is always larger than Mbar -except for 
IRAS 08339-1-6517, that has expelled a considerable fraction 
of its neutral gas to the i ntergalactic medium and shows a 
disturbed Hi kinematics (jCannon et al.ll2004 1 with a long 
tidal stream that makes impossible to get a good estimation 
of Mdyn (JLdpcz-Sanchcz, Estcban & Garcia-Roias 2006)-. 
Besides the uncertainties in Mdyn, this indicates the pres- 
ence of dark matter in all systems. The dark matter con- 
tribution would be even higher if, as we said, our values of 
Mdyn are underestimated because of the uncertainty in the 
extension of the H i disk (in all cases, except in those galax- 
ies for which interferometric data were available, we used 
the maximum of the radius of the optical extent to com- 



pute Mdyn)- The dotted yellow line in Fig (T^j indicates the 
position of Mdyn=Mbar if Mdyn is computed assuming that 
the extension of the neutral gas is 2.5 times the size of the 
optical extent. Indeed, only interferometric Hi maps and 
a detailed analysis of the rotation of the neutral gas (i.e., 
de Blok et al. 2008; Westmeier et al. 2010) can provide a 
better estimation of the dynamical masses of galaxies. This 
issue is even more important if interactions are disturbing 
the rotation pattern of the H i gas. A clear example of this 
is Tol 9 w ithin the Klemola 13 group. Our in terferometric 
H I map (jLopez-Sanchez et all l2008bl l2010b[ ) shows that 
the neutral gas cloud in which this BCG is embedded in- 
cludes not only Tol 9 but also some nearby dwarf galaxies. 
Indeed, this H i cloud seems to rotate as a single entity, 
and shows a long tidal tail in direction to other galaxies of 
the group. However, the maximum of the H i emission is 
located exactly at Tol 9. 

In any way, the important aspect to emphasize here is 
that, besides the unknown amount of dark matter, strongly 
interacting systems, as HCG 31 AC or Tol 1457-262, lie 
away from the observed main trend, having dynamical 
masses that are almost 14 and 7 times their baryonic 
masses. Mkn 1087, Mkn 1199 and Tol 9, that are in clear in- 
teraction with nearby objects, also have larger Mdyn than 
expected. Consequently, SBS 1054-1-365 and Mkn 5, that 
clearly have dynamical masses which are more than one 
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Fig. 10. Dynamical mass (Mdyn) vs. absolute B, V , R and J magnitudes for the galaxies analyzed in this work. A linear fit is 
shown with a dashed red line. 
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Fig. 12. Comparison between the baryonic mass (Mbar) and 
the dynamical mass (Mdyn) for our sample galaxies. The dashed 
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strongly interacting systems (HCG 31 AC or Tol 1457-262) lie 
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order of magnitude higher than those expected for dwarf 
galaxies with A/bar^ 10^ M0, may also have highly per- 
turbed H I kinematics. The same situation may be hap- 
pening in NGC 5253, which has a dynamical mass that is 
almost an order of magnitude higher than that expected 
for a galaxy with A/bar "^ 10^ Mq. 



6. Mass-metallicity relations 

Our data set allows to investigate the mass-metallicity 
{M — Z) relation of star-forming galaxies. The relation- 
ship between metallicity and stellar mass provides key 
clues about galaxy formation and evolution; however com- 
monly the luminosity is used instead of the mass to ana- 
lyze such correlations (i.e, Paper IV and references within). 
Observationally, the M — Z relation arises because low 
mass galaxies have larger gas fractions than higher mass 
galaxies (i.e., BoseUi et al. 2001; Kewley & Ellison 2008). 
Theoretically, the mean stellar metallicity of the galaxies 
increases with age as a consequence of the chemical en- 
richment of the ISM, while the stellar mass increases with 
time as galaxies undergo merging processes (i.e. Somerville 
& Primack 1999; Calura, Matteucci & Menci 2004). Once 
the NIR luminosities or the optical- iVXR SED are known, 
Mgtars can be estimated relatively well using stellar evolu- 
tionary synthesis models, as we explained in the previous 
section. Hence, the main problem to study the M — Z rela- 
tion lies in all the uncertainties involving the determination 
of an accurate oxygen abundance, such as different meth- 
ods yield to very different results (see Paper IV and Kewley 
& Ellison 2008). Here, the oxygen abundance of the major- 
ity of the galaxies was computed using the direct method, 
but Pilyugin (2001a) calibration has been applied to com- 
pute the metallicity of some few massive objects (Mkn 1087, 
Haro 15, IRAS 08339-1-6517, ESO 566-7, ESO 566-8), as we 
explained in Paper IV. 

Figure [T51 shows the relations between the stellar mass 
and the oxygen abundance, and Figure [TJ] shows the rela- 
tions between the baryonic mass (left panel) and the dy- 
namical mass (right panel) with the oxygen abundance. 
From Figures [T2] and [H] is quite evident that a. M — Z 
relation is satisfied for our sample galaxies. Although there 
is still a considerable dispersion for some objects, the com- 
parison with the luminosity-metallicity relation (see Fig. 17 
and Sect. 5 of Paper IV) suggests a closest correlation when 
using the stellar, baryonic, or the dynamical masses than 
the absolute optical/ A'XR magnitudes. 

The Mstars — Z diagram (Fig. [T31) shows a large 
dispersion for galaxies with 12-|-log(0/H)^8, as there 
are dwarf galaxies with Mstars ~ 4 x 10^ Mq 
(HCG 31 F and E, SBS 1054-^365) and large systems with 
Mstars - 5 X 10^ Mq (UM 420, Tol 1457-262 Obj 1) within 
this metallicity range. The origin of this dispersion is that 
the low-mass systems are TDG candidates, that have higher 
oxygen abundance than that expected for their mass (this 
is not the case of SBS 1054+365), while the high- mass ob- 
jects are very probably a merger of two independent galax- 
ies (and hence their oxygen abundance is much lower than 
the expected for a single, more massive galaxy). Neglecting 
the TDG candidates and the galaxies in the processing of 
merging, a linear fit to the data yields to 
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Fig. 13. Relation between Mstars and the oxygen abundance 
for our galaxy sample. A linear fit to the data is shown with 
a dashed red line. The TDG candidates are plotted with a red 
square, while galaxies in the process of merging are shown with 
a yellow diamond. A linear fit to the data without consider- 
ing these two groups of objects is shown with a continuous 
red line. Some previo us Mstars — Z relations a re also plotted: 
iTremonti et al.l J2004D with a yellow dotted line. lTremonti et al.l 
(|2004r ) corrected by a factor of 0.3 dex in oxygen abundance 
-dotted-dashed yellow l ine-, and the iKewlev fc EllisonI (|2008l ') 
relation considering the iPilvugiril (|2001br i empirical calibration 
to derive the metallicity -dashed pink line-. These relations are 
only valid for log Mstars > 8.5. 



being x=12-|-log(0/H). This relation has a correlation co- 
efficient of r=0.849 and it is plotted in Fig. [T3] with a 
red continuous line. Our Mstars — Z relation is quite dif- 
ferent from previous relations given in the literature. For 
comparison, the Mstars — Z d iagram includes the relations 
prov ided by Tremonti et all (|2004[ ) -do tted yellow line- 



and IKewlev fc EllisonI ()2008[ ) using the iPilvuginI ()2001b[) 
calibratio n to compute the oxy gen abundancqj -dashed 
pink line-. ITremonti et al.l (|2004 D derived the oxygen abun- 
dances using theoretical photoionization models, that over- 
estimate between 0.2 and 0.4 dex the metallicity derived 
from the direct T q method (see Paper I V). Hence, we also 
plot in Fig. [T^ the ITremonti et all ()2004|) relation corrected 
by a factor of 0.3 dex in oxygen abundance (dotted-dashed 
yellow line). As we see, ITremonti et al. (2004) relation is 
steeper than that derived here, and do es not agree well 
with our data. On the other hand, the IKewlev fc EllisonI 
(|2008l ) relation is quite flat in comparison with our obser- 
vational data. These authors commented that for masses 
between 3 x 10* and 10^^ Mq, the metallicity of their sam- 
ple galaxies rises only ^^0.2 dex on average, but we find a 
variation of ^0.7 dex in the s ame mass in terval. 

Interestingly, Kewley fc E llisonI (|2008f ) did not derive 
any Mstars ~ Z relation using oxygen abundances deter- 
mined with the To method, such as the SDSS catalog con- 
tains very few metal-poor and starbursting galaxies, and 
the scatter of the available data is huge. These authors 
finally concluded that the choice of the metallicity calibra- 



(6.03 ± 0.28) + (0.238 ± 0.031) logM^ 



^ As we concluded in Paper IV, this calibration provides the 
best results to the oxygen abundances derived for our galaxies, 



(31) that were mainly computed following the direct To method. 
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Fig. 14. Relation between A'har (left panel) and Mdyn (right panel) with the oxygen abundance for our galaxy sample. The TDG 
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tion has the strongest effect on the M — Z relation, be- 
catise a considerable variation in shapes and y-intercepts 
are found. Many of their fits suggested a flatter M — Z 
relation at hig her masses, some thing that was previously 
noticed by iTremonti et al.l (|2004f) . These authors explained 
this issue as a consequence of effective galactic winds that 
remove metals from the low-mass galaxies (M < lO^"'^ 
Mq). Although our data do not allow to explore this is- 
sue at high masses, we do not see any trend of this effect. 

It is very interesting to note that Tol 9, in which we 
detect a clear example of galactic wind, has a relatively 
low stellar mass for its expected metallicity. Probably that 
is indicating the strength and youth of the star-formation 
phenomena in this BCG. We should expect that the posi- 
tion of this object in the Mstars-O/H diagram will move to 
higher masses and lower metallicities if the star-formation 
processes continue and the fresh new material is expelled 
far from the galaxy via the effect of galactic winds. 

The linear fits to the Mt,ar — Z and M^yn—Z relations, 
that are plotted with red continuous lines in Figure [U and 
do not consider TDG candidates and mergers in progress, 
are 



X = (5.81 ± 0.25) + (0.375 ± 0.041) logMba^, 
X = (5.07 ± 0.31) + (0.291 ± 0.030) logMdyn, 



(32) 
(33) 



being x=12-|-log(0/H). The correlation coefficients are 
r = 0.821 and r = 0.929 for the Mbar - Z and Mdy„ - Z 
relations, respectively. The slopes of all the M — Z relations 
agree relatively well. We note the tightness of the M^yn—Z 
relation: except the mergers in progress (HCG 31 AC and 
Tol 1457-262), all galaxies are found relatively close to the 
relation. This indicates that the dark matter content of the 
galaxies also increases with the metallicity, in agreement 
with the predictions of the evolutionary galaxy models. 

In summary, the scatter observed in the luminosity- 
metallicity and in the mass-metallicity relationships of star- 
forming galaxies are consequence of both the nature and 
the star-formation histories experienced by these objects. 
Only a detailed analysis of each system can give the clues 
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Fig. 15. Relation between the SFR/area and the Hi gas den- 
sity for our galaxy sample. We show two values per galaxy , one 
assuming the Ha- based Sfr following the lKennicuttI (119981 ) cal- 
ibration (blue circles) and other considering the Sfr assumed 
combining all multi-wav elength data. The solid green line is the 
best fit to the M 51 data lKennicutt et aP (|2007l 1: the dashed red 
line is the relation for integrated values of star-forming galaxies 
derived by iKennicuttI (|l99^). 



needed to understand the evolution of the global proper- 
ties in star-forming galaxies and their comparison between 
dwarf, normal and massive galaxies. 



7. Schmidt-Kennicutt relation 

We now investigate if the studied galaxies do obey the 
Schmidt-Kennicutt scaling laws of star formation. It is well 
known that there exist a tight correlation between the av- 
erage Sfr per unit area and the mean surface density of 
the cold gas on global galactic scales. Such corr elation is 
usually pa,rameterized via a power-law relation ([Schmidtj 
ll959l . ll963l:[Kennicuttl[l998l:lKennicutt et al.ll2007ll that has 
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Table 8. Mass-to-light ratios of all mass estimations compiled in Table[7l The Mhai / Mdy-n, Mstars/Mbar, Mgas/Mtar, Mgas/Mstar 
and Mdust/Mgas ratios are also listed in the last columns. 



Galaxy 


Mhi 

Lb 


Mail 
Lb 


Afd„st 
Lb 




Mk^p 
Lb 


Mdyn 

Lb 


Lb 


Af.„. 
Lb 


Afdan 


Mt.r- 


Mga, 


Mg„e 


Mdust 


Mb.r 


M,t„„ 


Mgas 




r*f0i 


Tio* f 

L Lq 


1 TlO^fe 

J L Lq 


1 rio^fe] 

J L -^0 J 


r*foi 


r^oi 


\Me^ 


fAf©] 
UoJ 










[101 


HCG 31 


0.369 


3.86 


1.27 


18.7 






0.167 


0.655 




0.255 


0.744 


2.91 


2.60 


" AC 


0.395 


5.86 




19.8 


3.7 


9.23 


0.156 


0.678 


0.073 


0.231 


0.769 


3.34 




" B 


0.409 


0.96 




17.3 


0.55 


1.14 


0.169 


0.709 


0.622 


0.238 


0.762 


3.20 




" E 




7.09 




83.7 




4.15" 


0.150 














" F 


1.53 


13.7 




28.6 


0.75 


3.73 


0.097 


2.11 


0.565 


0.046 


0.954 


20.8 




" G 


0.350 


1.72 




12.6 


0.38 


1.23 


0.189 


0.651 


0.531 


0.290 


0.710 


2.44 




Mkn 1087 


0.138 


0.92 


0.69 


11.4 


0.50 


1.60 


0.256 


0.438 


0.274 


0.583 


0.417 


0.715 


3.78 


" N 




1.29 




12.8 


0.084 


2.78 


0.153 














Haro 15 


0.159 


1.79 


0.60 


22.2 


0.35 


1.05 


0.286 


0.495 


0.471 


0.577 


0.423 


0.732 


2.85 


Mkn 1199 


0.0419 


2.49 


1.07 


48.0 


0.28 


5.67 


0.945 


1.00 


0.177 


0.945 


0.055 


0.059 


19.2 


" NE 




0.78 




10.1 


0.17 


4.22° 


1.05 














Mkn 5 


0.273 


3.08 


3.76 


11.4 


8.0 


13.7 


0.575 


0.936 


0.068 


0.614 


0.386 


0.628 


10.4 


IRAS 08208+2816 




4.14 


1.73 


32.5 


0.076 


1.18" 


0.585 














IRAS 08339+6517 


0.0805 


2.19 


0.59 


34.2 


0.151 


0.56 


0.720 


0.827 


1.48 


0.871 


0.129 


0.148 


5.57 


" Comp 


0.234 


0.87 




1.37 




3.34 


0.623 


0.932 


0.279 


0.669 


0.331 


0.496 




POX 4 


0.219 


11.2 


0.183 


19.1 


0.098 


1.49 


0.281 


0.570 


0.383 


0.492 


0.508 


1.03 


0.63 


" Comp 




2.30 




36.4 






0.406 














UM420 




6.76 


4.9'' 


13.4 


0.204 


1.94° 


0.444 














SBS 0926+606 


0.664 


6.58 


1.39 


31.6 






0.285 


1.16 




0.246 


0.754 


3.97 


1.59 


" A 


0.746 


10.9 




28.0 




1.79 


0.226 


1.21 


0.676 


0.187 


0.813 


4.35 




" B 


0.587 


2.54 




34.5 




3.26 


0.340 


1.11 


0.342 


0.305 


0.695 


2.28 




SBS 0948+532 




31.7 




51.4 


0.57 


2.46° 


0.345 














SBS 1054+365 


0.929 


10.2 


1.83 


35.1 


1.19 


22.9 


0.505 


1.73 


0.076 


0.292 


0.798 


2.43 


1.50 


SBS 1211+540 


0.762 


6.61 




15.8 


3.57 


3.61 


0.321 


1.33 


0.368 


0.242 


0.758 


3.13 




SBS 1319+579 


0.408 


0.88 


0.75 


1.94 


2.14 


3.48 


0.669 


1.11 


0.318 


0.514 


0.486 


0.946 


1.40 


SBS 1415+437 


0.964 


7.35 




14.0 


2.50 


4.90 


0.465 


1.74 


0.355 


0.268 


0.732 


2.74 




III Zw 107 


0.378 


3.37 


0.71 


41.8 


0.046 


1.02 


0.349 


0.848 


0.834 


0.412 


0.588 


1.43 


1.42 


Tol9<= 


0.595 


4.32 


3.06 


34.8 


0.152 


7.38 


0.564 


1.35 


0.183 


0.418 


0.582 


1.39 


3.90 


Tol 1457-262 


0.258 


5.17 


1.01 


22.5 




5.22 


0.479 


0.820 


0.157 


0.584 


0.416 


0.711 


2.95 


" Obj 1 




5.70 




19.3 


0.051 


2.40° 


0.387 














" Obj 2 




4.25 




23.9 




2.49° 


0.601 














Arp 252 




3.01 


1.45 


18.9 






0.775 














" ESO 566-8 




3.57 




12.9 


0.21 


1.54° 


0.716 














" ESO 566-7 




1.09 




19.6 


0.047 


1.75° 


1.03 














NGC 5253 


0.070 


2.84 


0.182 


9.67 




3.60 


0.496 


0.594 


0.164 


0.841 


0.159 


0.189 


1.94 



Using the tentative value of M^y^ computed using Eg 1251 



proven to be very useful as an input scaling law for analyt- 
ical and numerical models of galaxy evolution (e.g., Kay et 
al. 2002). 

Figure [13] shows, on logarithm scale, the Sfr per unit 
area versus the surface density of the Hi gas (Mni/area) 
for all the galaxies for which we have Hi measurements. 
We plot two valu es for each g alaxy , one assuming the Ha- 
based Sfr using iKennicutn (|1998[ ) calibration (blue cir- 
cles) and other considering the Sfr assumed combining all 
multi- wavelength data (last column in Table [S]). Almost 
both values are quite similar in all galaxies except in some 
objects, remarkably POX 4 and NGC 5253. The major- 
ity of the galaxies are loc ated close to the relation given 
by iKennicutt et al.1 (|2007D . which is the best fit to the 
data of star-forming regions within the nearby Sbc galaxy 
M 51. These authors also included the molecular gas to get 
this relation, but we have not considered it in our galaxy 
sample. The assumption of neglecting the molecular gas 
is valid in low-mass, low-metallicity galaxies, because of 
both the difficulty of detecting CO and the uncertainties of 
the correspondence between CO and H2 in low-metallicity 
objects (i.e., Wilson 1995; Taylor, Kobulnicky & Skillman 
1998; Braine et al. 2004). However, we should expect some 



molecular gas contribution in more massive galaxies, as 
IRAS 08339+6517, Mkn 1087 and Mkn 1199. 

From Fig. [13 it is evident th at our data agree much 
better with the relation given by IKennicutt et al.l (I2007D 
than with the relation obtained bv iKennicuttI () 199811 for 
star- forming galaxies (and not regions within galaxies). 
Interestingly, a recent study of the star-formation activ- 
ity within UV-rich regions f ound in the outskirts of the 
galax y pair NGC 1512/1510 (jKoribalski fc Lopez-Sanched 
I2OO9D yield to this same result. We note t hat we plot the 
Ha-based S fr using t he 'Kennicutt! ("1998") calibrat i on be- 
cause both IKennicutt (1998) and Kennic utt et al.l (|2007f) 
relations use this calibration. I n case of usi ng the Ha-based 
Sfr derived from the .Calzetti et ahl (12007) calibratio n, the 
agreement of our data with the IKennicutt et aTJ ([2003) rela- 
tion will be even better. Some clear disagreements with the 
scaling laws of star formation are Mkn 1199 and NGC 5253 
(both seem to be H I deficient; the molecular gas component 
in Mkn 1199 would not explain its position in the diagram, 
as it would require that ^40% of the total neutral hydrogen 
mass is H2) and SBS 1319+579 (that show lower Sfr than 
that predicted considering their Hi gas amount). We then 
conclude that some external factors are indeed affecting the 
normal star-formation activity in these three galaxies. 



20 



Lopez-Sanchez: Massive star formation in Wolf-Rayet galaxies V: Multiwavelength global analysis 








iiii| 1 









^\^ HOG 31 F 










^\^ • 








1 


- SBS 1054+3b5>v# SBS 1415+437 


ff 




- 




• SB&-Bfl26+606A» 


S i 




- 




SBS 121 1+540 ^^^ 










SBS 0956*g06 B • 






- 




■ 


-!*• • 


9 •IIIZw107 


- 




MknS^ 


HCG31 
G31 G) 


»Tol 1457-262 


" 






m 


Haro15\»Mkn1087 




1 


_ 


A 


^\,,^^ 


_ 


















IRAS 08339+651 7 • ^\^ 






NGC5253^ ° 


^ 


-\^- 






^ 




"^ 








• Mkn1199 


- 




Linear fit: y = ( 2.32 + 0.41 ) + ( 


0.309 + 0.045 
.,..1 


)x. r=0.858 
1 1 





1E7 1E8 1E9 1E10 



1E11 







1 1 1 1 

Linear fit: y =( -2.51 ± 0.38 ) - ( 0.01 6 ± 0.041 )x, r=0.0683 




■ 


0.01 


- 


• HCG 31 E 






- 






SBS 1054+365 

• • 
• 


SBS 0948+532 

• III Zw 107 

■ , Tol9,* 

Hnn 31 AC • «. 


Mknn99 

* IRAS 
# •08339+6517 
IRAS 08208+2816 


" 




SBS 1415+437 


. • •• * 


•Arp 252 


1E-3 


r 


SBS 121 1+540 • 

• 
Mkn5 


• • • 

* ^ UM 420 
NGC52?3 Mkn1199NE 


•Mkn 1087 


- 








SBS 1319+579^ 

^ IRAS 08339+6517 


comp 


" 


1E-4 






1 




" 



1E7 1E8 1E9 1E10 1E11 

M, [Me ] 

Stars ■■ ^ ■■ 



3E-5 



■ 





'"! 1 


' ■ ■ 






♦ SBS 0948+532 


■ 




"™^"=» SBS 0926+606 A 
SBS1054+365» • 


• POX 4 






---._._^ » SBS 1415+437 
SBS 1211+540"^"" — -~~^_^ 


UM 420 
• • 

• • Tol 1457-262 
HCG 31 AC • 




Mkn5 -- ^_ ^•mZwIOT 

POX 4 Comp • SBS 0926*^606 B ^°^ '''' ----_M1''__1 

• HCG31G ,,*., 
Haro15 


#IRAS 08208+2816 . 
• 

• Arp 252 

"~~~-^&J36339+6517- 




Mkn1087N» 






T 


HCG 31 


C SBS 1319+579 


•Mkn 1087 ~ 




Mkn 1199NE 


: 




Linear fit: y = ( -2.20 + 0.54 ) + (-0.1 40 + 0.059)x. r=0.388 

1 1 1 1 


■ 



1E8 1E9 1E10 

M, [M„] 

Stars •■ ©J 



0.2 





' ' ' ' 


■ 




• SBS 1054+365 

• Mkn 5 


- 


- 


HOG 31 AC 

♦ Tol 9 


- 




■ 


" 




SBS 1415+437 Tol 1457-262 ^ 

• ♦ ^ Mkn 1199 


■ 




SBS 1211+540 NGC5253 SBS 1319+579 
HCG 31 F ggg 0926+606 B IRAS 08339+6517 Comp 


■ 




SBS 0926+606 A • 

»Dnv^ » Mkn 087 
HCG 31 G« POX 4 


" 


r 


"^^^^ III Zw 107 


- 




■ 


■ 




IRAS 08339+6517 






♦ Merger in progress 






■ Very perturbed Ha and/or H 1 kinematics 





1E7 1E8 1E9 1E10 

M.„„. [Me] 



1E11 



Fig. 16. Comparison between the stellar mass and some mass-to-light ratios for our sample galaxies. 



8. Analysis of the mass-to-light ratios 

Table \E\ compiles all the mass-to-B luminosity ratios de- 
rived in this work. Some interesting relations are plotted 
in Figure [ini that compares some mass-to-light ratios with 
the stellar mass derived from the i7-band luminosity. 

The H I mass-to-light ratio of a galaxy is a distance- 
independent quantity that compares the H I mass with 
the luminosity in the B-band. This property correlates 
with many galaxy par ameter, as the galaxy typ e, galaxy 
color or galaxy mass (jRoberts fc Havnesl |1994[ ). Indeed, 
the comparison of the Mm/ Lb ratio with the stellar 
mass in our sample galaxy clearly indicates that less 
massive galaxies have a higher mass fraction of neutral 
gas. The majority of the galaxies have a H i-mass-to- light 
ratio between 0.1 and 1.0, in agreement w ith previous 
estimations in star-forming dwarf galaxies (Salzer et al.l 
I2OO2I: iHuchtmeier. Krishna fc Petrosianl 120051 ). We note 
some peculiar objects in this diagram. The H i-mass of 
Tol 9 has been overestimated because the H i cloud in 
which it is embedded includes s everal dwarf galaxies 
([LoDez-Sanchez et all IJoOSbl l2010b[ ). On the other hand, 
two galaxies (Mkn 5 and NGC 5253) are very Hl-deficient. 
In particular, NGC 5253 is very far from the typical posi- 
tion of the galaxies, showing a M-hi/Lb of ^0.051 Mq/Lb- 
The M}ii/Lb ratio of Mkn 1199 is also slightly low, even 
for a massive galaxy. As we already suggested (Sect. 3.4.3 of 
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Fig. 17. Comparison of the Mgas/Mstars ratio with the oxygen 
abundance for our sample galaxies. The dotted horizontal pink 
line indicates the position of Mgas=M stars- The red continuous 
line is a linear fit to the data, neglecting Tol 9, for which the 
neutral gas mass corresponds to this object and some dwarf 
surrounding galaxies (see text and LS08). 
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Paper II), Mkn 1199 may has lost part of its neutral gas in 
the interaction process with its NE companion. Neglecting 
the contribution of Tol 9, Mkn 5 and NGC 5253, a linear 
fit provides the empirical relation 



log 



Mhi 



= (2.32 ± 0.41) - (0.309 ± 0.045) log A4 



(34) 



that has a correlation coefficient of r = 0.858. 

We do not find any high M-hi/Lb ratio (> IMq/Lb) 
in our sample galaxy, except in the case of the TDG 
candidate HCG 31 F, that has 1.53 Mq/Lb- High Hi 
mass-to-light ratios has been reported in some few galax- 
ies. The detailed analysis of the gas-rich low surface 
brightness dwarf irregular galaxy ESQ 215 -G009 per- 
formed by IWarren. Jerien &: Koribalskil ()2004l) confirmed 
an extremely high Mm/ Lb of 22 ± 4 Mq/Lb in this 
galaxy, for which the H i disk extends 6.4 ± 0.4 times 
the Holmberg radius. They concluded that ESO 215-G009, 
that is very isolated (no neighbors identified out to 1 
Mpc), has a low Sfr that probably remained unchanged 
thro ughout the galaxys existence. In a subsequent pa- 
per (jWarren. Jerien fc Ko ribalski"2006') these authors sug- 
gested that high Mm/ Lb galaxies are not lacking the 
baryons to create stars, but are under luminous as they lack 
either the internal or external stimulation for more exten- 
siv e star formation. 

IWarren. Jerien fc Koribalskil (|2003) derived an empiri- 
cal upper envelope for Mm/ Lb as a function of the abso- 
lute B-magnitude, that accounts for the maximum amount 
of atomic hydrogen gas a galaxy of a particular luminosity 
can retain in the Universe today. All our sample galaxies 
satisfy this empirical relation. 

The derived Myhi/Lb ratios for our galaxy sample lie 
in the range lO^'^-lO^'' Mq/Lb- Although the scatter of 
our data is high (a tentative linear fit gives a very low cor- 
relation coefficient), they indicate that the Mhii/Lb ratio 
slightly decreases with increasing stellar mass, suggesting 
that the ionized gas to stars ratio is higher in dwarf galax- 
ies. SBS 0948-1-532 is away from the rest of the objects 
because of its very high Ha fiux (see Sect. 3.11 in Paper I). 
On the other hand, the ionizing star cluster mass-to-light 
ratio. Mi,/ Lb, seems to be rather constant with the stel- 
lar mass, showing an average value of '^0.0022 Mq/Lb- 



Two galaxies, SBS 1319-1-579 and the companion object of 
IRAS 08339-1-6517 lie apart from this tendency, as we also 
saw in Fig. [9] 

The M^yn/LB ratio seems to slightly decrease with the 
stellar mass. However, the analysis of this diagram is dif- 
ficult, because interactions notably modify the estimation 
of the dynamical mass. Usually, perturbed kinematics yield 
to higher Mdyn (HGC 31 AC, Mkn 1199, Tol 1457-262, 
Tol 9), but sometimes the existence of tidal tails with a 
rather constant velocity give a lower Mdyn than the real 
one (IRAS 08339+6517). SBS 1319+579 (a probable merg- 
ing of two dwarf objects) and IRAS 08339+6517 Comp (in 
interaction with the main galaxy of the system) also show 
somewhat high Mdyn/is ratios. Hence, we may suggest 
that galaxies Mkn 5 and SBS 1054+365, that lie far from 
the rest of the objects, have a perturbed kinematics, being 
-Mdyn overestimated in both cases. 

Figure [T7| compares the ratio between the gas and 
the stellar masses with the oxygen abundance. Clearly, 
M gas /M stars dccrcascs with increasing mctallicity, indicat- 
ing that the importance of the stellar component to the 
total mass is higher in more massive galaxies. We should 
expect this result, which is qualitatively the inverse behav- 
ior of the M-hi/Lb ratio with the stellar mass. A tentative 
linear fit to the data (excluding Tol 9) is plotted in Fig.[T7| 
with a continuous red line, and provides 



log 



A^a 



gas 



Ms 



(11.7 ±3.2)- (1.42 ±0.39)a;. 



(35) 



being x=12+log(0/H), and has a correlation coefficient of 
r — 0.655. Following this analysis, we should expect that 
galaxies with 12+log(0/H)r-^8.2-8.3 have relatively equal 
gas and stellar masses. HCG 31 F is the galaxy with highest 
M gas /M stars ratio (and lowest M stars/ Lb ratio), indicating 
its low content of evolved stars. This result agrees with 
our suggestion (Lopez-Sanchez ct al. 2004a) that this TDG 
has created stars mainly using the ne utral gas from the 
long arm-like H i structure (|Verdes-Montenegro et al.ll2005() 
found between members AC and G of the HCG 31 galaxy 
system. 

Finally, we compar e some mass-to - light ratios with the 
colours of the galaxies. lAmorin et al.l ()2009[ ) reported that 
the underlying component (host) of blue compact galax- 
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ies is redder with decreasing M-hi/Lb- We do not find 
any correlation between the Mm/ Lb ratio and the opti- 
cal colors of the underlying component, but we did not 
perform a detailed analysis of the str uctural parameters o f 
the host underlying the starburst as lAmorin et alj ()2009[ ) 
did. The comparison of the global B — R colour and the 
neutral gas mass-to-light ratio is shown in the left panel of 
Fig HI and has a huge scatter. We should remember that 
Mkn 1199, Mkn 5 and NGC 5253 seem to be Hi deficient, 
and th at IRAS 08339-1-6517 is a luminous blue compact 
galaxy (jLopez-Sanchez et al.ll20Q6[ ). and hence their real po- 
sitions in this diagram is uncertain. Although a tentative 
fit to the data suggests that galaxies with redder B — V 
colours have lower Mm/ Lb ratios, the huge scatter does 
not allow us to confirm such tendency. 

However, we do observe a clear relation between the 
stellar-to-light ratio and the global B ~ R colour of the 
galaxies (right panel of Fig. [TS|) . This tendency seems 
to be also a consequence of the building of the galax- 
ies, as more massive galaxies have experienced more star- 
formation events than less massive objects, and hence tend 
to show redder stellar populations than dwarf galaxies. This 
result also agrees quite well with the observed tendencies 
that Mm/ Lb ratio decreases with the stellar mass and that 
Mgas/Mstars ratio decreases with increasing metallicity. 

9. Dust properties in star-forming galaxies 

Our data set allows us to investigate the properties and 
effects of the dust content in low-metallicity star-forming 
galaxies. Figure [12] plots the reddening coefficient, c(H/3) 
-obtained using our optical spectra-, as a function of 
the warm dust mass, Afdust -derived from FIR data-. 
Neglecting the data for UM 420 (as we said before, Afdust 
is overestimated because of the FIR contribution of the 
foreground galaxy UGC 01809), we see a clear correlation 
between both quantities: galaxies with higher amount of 
warm dust (and hence, as we saw in Figure |9l higher lumi- 
nosity) show higher extinction. This conclusion is in agree- 
ment with other results previously found in this work, as 
the correlation between c(H/3) and the oxygen abundance 
discussed in Paper IV. More important, this result indi- 
cates that most of the dust is internal to the galaxy and 
not in the line of sight. Detailed analysis of the dust dis- 
tribution within nearby galaxies (i.e. Muiioz-Mateos et al. 
2009) found clear relationships between the dust content 
and general properties of nearby spiral galaxies, such as 
galaxy type, luminosity, and metallicity. Here we confirm 
in an independent way, such as the extinction was derived 
from our optical spectra, that the dust content and there- 
fore the extinction in dwarf galaxies depends on their metal- 
licities and luminosities, and very probably also on their 
star-formation histories. A proper estimation of the amount 
of dust within such objects is needed to perform appropri- 
ate statistical analysis involving larger galaxy samples. 

We now investigate the dust-to-gas ratio, Mdust/Mgas, 
of our sample galaxies, a very important quantity when 
studying the chemical enrichment of the ISM, as it accounts 
for the amount of metals locked up onto dust grains through 
the stellar yields. The existence of a correlation between the 
Mdust/Mgas and the oxygen abundance has been reported 
in many studies (i.e., Lisenfeld & Ferrara 1998; James et 
al. 2002; Draine et al. 2007; Muiioz-Mateos et al. 2009). 
Figure ^U\ shows the dust-to-gas ratio as a function of the 
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Fig. 19. Reddening coefficient, c(H/3), vs. warm dust mass, 
Mduat. We note that the Mdust value in UM 420 is overestimated 
because of the FIR contribution of UGC 01809. 



oxygen abundance. The gas mass was computed assuming 
only the Hi and the Hei gas, but not the contribution of 
the molecular gas, that as we already explained is not im- 
portant in dwarf low-metallicity objects. The derived dust- 
to-gas ratio of each galaxy for which we have both H i and 
FIR data are compiled in the last column of Table \E\ From 
Figure [20l it is evident that objects with higher metallic- 
ities tend to have higher Mdust/Mgas ratios, such as the 
amount of dust increases while the neutral gas is consumed 
as the galaxies experience new star-formation phenomena. 
The linear fit to our data (continuous red line in Fig. [201) 
provides this tentative relation 

\og{Mdust/Mgas) = (-12.0 ± 2.9) + (1.02 ± 0.36)a;, (36) 

with x=12+log(0/H) and correlation coefficient of r = 
0.637. Mkn 5 lies far the the majority of the points, but 
as we already commented this object seem to be very de- 
ficient in H I, so we did not include this point in the fit. 
Although with higher uncertainties, we also observe that 
more massive objects also tend to have higher Mdust/Mgas 

ratios. 

iDraine et al.l ()2007[ ) provided a relation between 
Mdust/ Mgas and the oxygen, abun dance -computed follow- 
ing the iPilvugin fc ThuanI (|2005[ ) calibration- in a sample 
of spiral and irregular galaxies, that we may re-write as 



log{Mdust/Mgas) = GAS + X, 



(37) 



being x=:12+log(0/H), and plotted in Figure [501 with 
a dotted pink line. The factor 6.48 was derived 
from \og[{Mdust/Mgas)MW /'i^-52] + xmw assuming 
[M dust /M gas) MW = 0.010 and xmw = 8.6 for the Milky 
Way. As we see, this relation lies away fr om our data, 
although the slope is the same in both cases. iDraine et al.l 
(2002) also found that the global dust-to-gas ratio of 
all their galaxies with 12-|-log(0/H)<8.1 falls below 
this equation, sometimes by a factor larger than 10. 
The dashed-dotted green line in Figure [50] plots the 
[Draine et al.[ (|2007[ ) relation divided by a factor 10. As 
we can see, t his re lation agrees better with our data. 
IDraine et al.[ ()2007l) also remarked that many of the 
low-metallicity galaxies have large H i envelopes mainly 
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Fig. 20. Dust-to-gas ratio, Mdust/Mgas, vs. the oxygen abun- 
dance for our sample galaxies. The red continuous line is a fit 
to our data (neglecting M kn 5). The dashed yellow l ine indi- 
cates the relation found bv lMufioz-Mateos et alj ([2009) analyz- 
ing the radial dust-to-gas profiles for a larger sample of spiral 
galaxies. This relation has been corrected by 0.3 d ex be cause 
it was derived assuming the iKobulnickv &: KewlevI (|2004l ) cali- 
bration to compute the oxygen abundances, that overestimates 
in 0.2-0.4 dex the oxygen abundance provided by the direct 
Te method. The dotted pink line is the relation provided by 
iDraine et al.l (2007) in their analysis of a sample of spiral and 
irregular galaxies. The dashed-dotted green line indicates the 
IDraine et al.l (|2007l ) relation divided by a factor 10. 



composed by un-enriched material. However, the metallic- 
ity is derived from the brightest H ii regions, that usually 
show the highest metallicities (most recent star formation) 
within the system. Consequently, the derived Mdust/Mgas 
ratio should correspond to lower oxygen abundances. 

Recently, iMunoz-Mateos et al.l (|2009D analyzed the ra- 
dial dust-to-gas profiles for a larger sample of spiral galax- 
ies and found a steeper relation between Mdust/Mgas and 
the metallicity, that they explained because the outskirts 
of spiral galaxies seem to have a much lower Mdust/Mgas 
than the central regions (see their figure 16). These au- 
thors suggested a link with the behavior found in dwarf 
galaxies because in the external regions of the spiral galax- 
ies the neutral gas has not yet undergone star-formation 
activity. However, this trend could also be a consequence of 
the radial decrease of the star-formation efficiency found in 
nearby spirals (i.e., Thilker et al. 2007; Leroy et al. 2008). 

We may check these hyp otheses comparing our dat a 
with the relation provided by iMufioz-Mateos et al.l ()2009D . 
However, we should modify slightly their equation, be- 
cause these authors used the IKobulnickv fc KewlevI (|2004D 
method to derive the metallicities, and this calibration over- 
estimates the oxygen abundances provided by the direct T^ 
method in 0.2-0.4 dex (see Paper IV). The modified rela- 
tion that is plotted with a dashed yellow line in Figure [^ 



\0g{M dust /M gas 



-23.77 + 2.45(2: -f 0.3), 



(38) 



being x=12-flog(0/H). We see that except for Mkn 5 all 
our data have lower Mdust/Mgas ratios that those pre- 
dicted following this relation. Hence, we suggest that the 
low Mdust/Mgas ratios found in dwarf low- metallicity galax- 
ies is a consequence of the large reserves of un-enriched 



neutral gas in the se systems, while the low M dust/Mgas 
ratios reported by 'Muno z-Mateos et al.l (120091 ) in the ex- 
ternal regions of spiral galaxies are due to the decreasing 
of the star-formation efficiency in those areas. Indeed, the 
outskirts of spiral galaxies seem to show lower H i surface 
densities -log(MH /area) ~0.4-0.8 Mp, p c~^ for external 
regions of M 51 (^Kennicutt ct al.' l2007li and NGC 1512 
(Koribalski fc Lop cz-Sanchcz 2009)- than the neutral gas 
envelopes of dwarf galaxies -log(MH/area) '--^l. 4-2.0 Mq 
pc~^ following Fig. [T51 -. that is translated in a decreas- 
ing of the star-formation activity following the Schmidt- 
Kennicutt law. However, more high-quality data and a de- 
tailed analysis of the dust and the H i distribution in nearby 
dwarf galaxies should be performed to confirm all these is- 
sues. 



10. Comparison with the closed-box model 

To study the environment effects of the gas content and 
the chemical enrinchment in galaxies, it is common to com- 
pare with the so - called closed-box chemical evolution mode l 
(iSchmidtl [1961 iSearle fc Sarge^ [1971 lEdmundsl [T990l) . 
According to this model, a galaxy consists initially of gas 
with no stars and no metals. The stellar IMF is assumed to 
be constant on time. Stars that end their life as supernovae 
are assumed to enrich the ISM with metals immediately. 
Throughout its life, the galaxy experiences instantaneous 
recycling and the products of stellar nucleosynthesis are 
neither diluted by infalling pristine gas nor lost via outfiow 
of enriched gas. Hence, the metallicity at any given time is 
only determined by the fraction of baryons which remains 
in gaseous form. The model can be written as 



Zo = yoln(l/Ai), 



(39) 



where Zq is the oxygen mass fraction, j/o is the yield by 
mass and /z is the ratio of the gas mass to the baryonic mass, 
/i = Mgas/Mbar- The gas mass corresponds to the hydrogen 
atomic gas with a correction for neutral helium, but it does 
not include molecular gas (Mgas = 1.32A/Hi) that, as we al- 
ready said, can be neglected in low-metallicity galaxies. 

Left panel of Figure [Hj compares the observed oxy- 
gen abundances to those predicted by closed-box mod- 
els, that are plotted with lines with different j/o- The 
green continuous line indicates the model with j/o=0.0074, 
that is the theoretical yield of oxygen expected for stars 
with rotation following Mevnet fc Maeden ([20021 ) models 
([van Zee fc Havnea 120 06^ . As we can see, the majority of 
the galaxies show oxygen abundances lower than the ex- 
pected by the closed-box models. The yield of oxygen that 
best fits our data (the effective yield) is yo=0. 003-0. 005, 
in agreement with previous results found in the literature 
(i.e., Lee et al. 2003; van Zee fc Haynes 06; Lee, Zucker fc 
Grebel 2007). Hence, the sample galaxies are generally not 
well reproduced by the simple closed box model, and there- 
fore inflow of pristine gas or outflow of enriched gas have 
played an important role in their chemical evolution. 

Interestingly, there is a galaxy that appears far away 
from the predictions given by the theoretical closed box 
model, but in the opposite direction to the rest of the 
galaxies. Indeed, HCG 31 F shows a much higher oxy- 
gen abundance than that expected following the closed 
box model. The explanation of this behavior is that this 
object is a TDG that was very probably formed from 
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Fig. 21. (Left panel) Comparison of the observed oxygen abundance and that predicted by simple closed-box chemical evolution 
models with instantaneous recycling and constant star-formation rates. The continuous g r een line indicates the ex pected trend 
if the galaxies are closed boxes with an oxygen yield t/o=0.0074 (jMevnet fc Maeden 120021 : Ivan Zee fc Havnesll200q ). Closed-box 
models with j/o =0.0015, 0.003, 0.005, and 0.03 are also plotted. A linear fit to our data is shown with a continuous red line, 
and indicates an oxygen yield of j/o =0.003-0.005. (Right panel) Effective yield plotted as a function of the Mgas/ Lb ratio. The 
closed-box yield is plotted by a dotted yellow line. A linear fit to the data is shown with a continuous red line. 



the material stripped from HCG 31 AC during the fly- 
by encounter between member G and the A-l-C com- 
plex (Lopez-Sanchez et al...2004a) . The TDG has acretted 
a large fraction of the pre-enriched H i gas available in 
the arm-like structure and now hosts a very intense star- 
formation activity. 

Finally, the comparison of the effective yield derived in 
each object with some global galaxy parameters (dynam- 
ical and baryonic mass, absolute magnitude, gas mass-to- 
luminosity ratio, and surface star formation rate) does not 
show any clear tr end. This result is almos t the same than 
that observed by Ivan Zee fc Havnea (|2006ll in the analysis 
of a sample of isolated dwarf irregular galaxies. The dif- 
ference is that these authors reported a strong correlation 
with the gas mass-to-luminosity ratio, that they explained 
in the sense that gas-rich galaxies are more likely to be 
closed boxes. But, as we see in the right panel of Figure ETl 
such tight correlation is not satisfied by our data, although 
we do observe the trend that galaxies with higher Mgas/ Lb 
ratios have lower effective yields (a linear fit to our data is 
shown with a continuous red line). Therefore, for intense 
star-forming and gas-rich galaxies, the closed box model is 
also not valid. We then conclude that environment effects 
are playing a crucial role in the evolution of these galaxies. 

11. Quantification of tlie interaction features 

Throughout this paper series we have compiled new ev- 
idences of the interaction-induced star-formation activ- 
ity in starburst galaxies, in particular in dwarf galax- 
ies. Alternative mechanism s, as the Stochasti c Self- 
Prop agating Star Formation (jGerola. S eiden. & Sc hulmanI 
[1983) model (that assumes statistical fluctuations of Sfr), 
or the ideas of the cyclic gas re-processing of the ISM 
(|Dave fc Oppenheimei|[l988l ) or gas compression by shocks 
due to the mass lost by galactic winds followed by the cool- 
ing of the ISM (Thuan 1991; Hirashita 2000 ), fail to explain 
some observational characteristics and the triggering mech- 
anism of dwarf starburst galaxies (i.e., Garcia-Lorenzo et 



al. 2008; Cairos et al. 2009). In the previous section, we 
demonstrated that the closed-box model is not valid to ex- 
plain the chemical evolution experienced by our sample of 
galaxies, emphasizing the idea that environment effects are 
needed to understand their observed properties. Indeed, the 
interaction/merger scenario explains, in a natural way, the 
starburst activity in these objects just as a consequence of 
the evolution of the galaxies throughout the cos mic time fol- 
lowing hierarchical formation models (IKa uffmann fc Whitd 
119931 ; (Kaufmann et all 119971; ISpringel et al.i [200^. These 
models predict that most galaxies have formed by merg- 
ing of small clouds of protogalactic gas and that galaxy 
interactions between dwarf objects were very common at 
high redshifts. 

However, the interaction features in dwarf ob- 
jects are, in many cases, not evident because of 
the lack of deep and h igh-resolution images and 
spectra (jMendez fc EstebanI [2000) and detailed multi- 
wavelength analyses. Now it is well known that inter- 
actions in dwarf galaxies are not usual with nearby 
giant galaxies (ICa mpos-Agu ilar. Moles fc Maseeosal 119931 : 
iTelles fc Terlevichill995t iTelles fc M addox"2000') but with 
low s u rface brightness galaxies ([Wilcots. Lehman fc Millei] 
Il99fil : iNoeske et all 120011: iPustilnik et all I2001D . or Hi 
clouds (e.g., Taylor et al. 1993, 1995, 1996; Thuan et 
al. 1999; van Zee, Salzer fc Skillman 2001; Begum et al. 
2006; Ekta, Chengalur fc Pustilnik 2006; H utchmeier et al. 
2008; Lopez-Sanchez fc Esteban 2008). Mendez fc Esteb"aiil 
(|2000D suggested, for the first time, that interactions with 
or between dwarf objects could be the main star for- 
mation triggering mechanism in dwarf galaxies. Later, 
lOstlin et all (|200l[ ) and lBergvall fc OstlinI (|2002D suggested 
that a merger between two galaxies with different metal- 
licities or infall of intergalactic clouds could very proba- 
bly explain the starburst activity in the most luminous 
BCDGs. Since then, studies focused on individual ob- 
jects have also shown that interactions d o play a deci- 
sive role in the evolution of these systems ([Johnson et al.l 
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Table 9. Interaction features in our WR galaxy sample. 
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" Some interaction features in this galaxy were previously reported by other authors. 

^ The chemical differences a nd the kinematics features can be explained by other reasons. See lLopez-Sanchez et aU | |2007H and 
iLopez-Sanchez et all ll2010af) . 
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Our exhaustive multi-wavelength analysis of starburst 
galaxies combining broad-band optical/ ATi? and Ha pho- 
tometry, optical spectroscopy, and X-ray, UV , FIR, 21-cm 
Hi line, and 1.4 GHz radio-continuum data compiled from 
the literature allows us to perform a quantitative analysis 
of the interaction features detected in each object. A sum- 
mary of the results found in each individual system of our 
WR galaxy sample is presented in Appendix A. To quan- 
tify the interaction features, we compile in Table [3] some 
interaction indicators classified in several categories, that 
we describe below. 

1. Morphological features, as the detection of faint 
plumes or bridges (Haro 15, IRAS 08339+6517, 
SBS 0926+606 B, SBS 1211+540, HI Zw 107, Tol 9, 
Tol 1457-262, Arp 252), prominent tails (HGG 31, 
Mkn 1087, IRAS 08208-1-2816, UM 420, SBS 0948-1-532, 
Arp 252), disturbed morphology (HGG 31, POX 4, 
Tol 1457-262), TDGs candidates (HGG 31, Mkn 1087, 
SBS 0926-^606 B, Arp 252) or mergers. 

2. Kinematical features detected in the analysis of the ion- 
ized gas (see Paper II) and the neutral gas (only in those 
systems for which interferometer H i maps are avail- 
able). The kinematical evidences found in the ionized 
gas of our sample galaxy includes: presence of objects 
with velocities decoupled from the main rotation pat- 
tern (Mkn 1087, Haro 15), sinusoidal velocity patterns 
that suggest a merging process (HCG 31 AG, Mkn 1199, 
IRAS 08208+2816, SBS 0926-1-606 A, HI Zw 107, 
Object 1 in Tol 1457-262), reversals in the velocity dis- 
tribution (Tol 9, Arp 252), indications of tidal streaming 
(HGG 31, IRAS 08208+2816, SBS 1319+579, Tol 9), or 
the presence of TDG candidates (HGG 31 Fl and F2, 
Mkn 1087, IRAS 08339+6517, POX 4, Tol 1457-262). 



3. Ghemical abundance differences within several star- 
forming regions within the same system: Mkn 1087, 
Haro 15, and Mkn 1199 are clearly interacting with 
dwarf galaxies with lower 0/H and N/0 ratios. 
NGG 5253, IRAS 08208+2816, and Tol 1457-262 con- 
tain zones of different chemical composition. In the case 
of NGG 5253, this is produced by localized pollution of 
massive stars, but in the cases of IRAS 08208+2816 and 
Tol 1457-262 the different chemical compositions seem 
to be caused by the regions corresponding to different 
galaxies in interaction. 

4. Furthermore, our multiwavelength analysis has pro- 
vided us further indications of galaxies that do not 
follow their expected behavior. The analysis of the 
mass-to-light ratios indicates: very low Mm/ Lb in 
SBS 1319+579 and NGG 5253, high Mhh/Xs in 
SBS 0948+532, low IVh/Le in SBS 1319+579, and 
high Mdyn/Ls in HGG 31 AG, Mkn 1199, Mkn 5, 
SBS 1054+365, Tol 9, and Tol 1457-262. Other ev- 
idences are: low H i mass content from single-dish 
data (Mkn 1199, Mkn 5), very extended Hi emis- 
sion embedding several nearby galaxies as HGG 31 
(VM03) and Tol 9 (LS08, LS+lOb), high Mxcp (Mkn 5, 
IRAS 08339+65 Gomp, SBS 1211+540, SBS 1319+579), 
or important deviations of the star-formation law 
(Mkn 1199, SBS 1319+579 and NGC 5253). These fea- 
tures may have been produced by interactions (loss of 
H I mass, enhancing of the star-formation activity, per- 
turbed dynamics) but they are just indirect evidences 
that should be confirmed by new deep observations (i.e.. 
Hi maps). 

A question mark in Table [3] indicates that the available 
data do not allow us to confirm this indicator. Last col- 
umn in Table [9] compiles the interaction degree that each 
system is experiencing after considering all positive inter- 
action indicators. We have divided the interaction degree 
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in four classes: low (no clear signs of interactions), proba- 
ble (there are some interaction indicators, but deepest data 
are needed to confirm it), high (we found clear evidences of 
interactions, but we do not see merger features) and very 
high (in the cases of finding clear merger features). 

Evident mergers between independent objects with rel- 
atively similar masses (major mergers) have been detected 
in HCG 31 AC, SBS 0926-1-606 A, IRAS 08208-f2816 and 
Tol 1457-1-262. The galaxy pair Arp 252, that is composed 
by ESO 566-7 and ESO 566-8, seems also to be experienc- 
ing the first stages of a major merger. Minor mergers are 
found in Haro 15 and Mkn 1199. UM 420 seems also to 
be experiencing a merger, as deep 3D optical spectroscopy 
(James et al. 2009) suggests. All these galaxies have a very 
high interaction degree. 

Mkn 1087, IRAS 08339-^6517, POX 4, SBS 0926+606 B, 
III Zw 107 and Tol 9 are experiencing clear interactions 
with nearby dwarf objects. In the case of POX 4, we still 
have to investigate ()L6pez-Sanchez et al.ll2010bl ) if its dwarf 
companion galaxy is a TDG candidate and the interaction 
was with a nearby diffuse H i cloud, or if this object actu- 
ally is an indep endent dwarf galaxy that crossed the main 
body of POX 4 (jMendez fc Estebanl[T997h . These six galax- 
ies have a high interaction degree. 

On the other hand, we find probable evidences of in- 
teraction in SBS 0948-1-532 (enhanced star formation ac- 
tivity, long optical tail, probable disturbed kinematics of 
the ionized gas), SBS 1211-t-540 (diffuse optical plumes, mi- 
nor merger indications), SBS 1319+579 (peculiar Mhi/Lb, 
Mmi/ Lb, and M-^/Lb ratios, perturbed kinematics sug- 
gesting merging or tidal stream phenomena). The Hi data 
of NGC 5253 shows disturbe d morphology and kinematics 
(jLopez-Sanchez et al.ll2008al ). suggesting that this BCDG 
has dis rupted or accreted recently a dwarf gas-rich com- 
panion (jKobulnickv fc Skillmanll2008l : iLopez-Sanchez et al.l 

l2ninaD . 

Only Mkn 5, SBS 1054+364, and SBS 1415+437 do 
not show interaction evidences in our exhaustive multi- 
wavelength study. However, Mkn 5 seems to be H i-deficicnt 
and seems to possess a perturbed neutral gas kinematics be- 
cause of its relatively high Mdyn/is ratio. SBS 1054+364 
also show a high Mdyn/^s ratio that may suggest per- 
turbed H I kinematics. Furthermore, the chemical abun- 
dance of this galaxy is very high in comparison to that 
expected from its baryonic mass. On the other hand, 
SBS 1415+437 shows a relatively low oxygen abundance 
for its baryonic mass. 

Considering all indicators, we find that 13 up to 20 
systems (68% of our WR galaxy sample) are classified 
with a high or very high interaction degree. We note, 
however, that four of these objects (HCG 31, Mkn 1199, 
IRAS 08339+6517 and Arp 252) show weU known inter- 
action evidences, but our analysis reinforces the evidences 
and improve our knowledge of these systems. Only three 
galaxies (Mkn 5, SBS 1054+364, and SBS 1415+437) do 
not show interaction features, but they show consider- 
able divergences of some properties when comparing with 
similar objects. Hence, it is evident that the majority 
of the analyzed galaxies (17 up to 20) are interact- 
ing or merging with or between dwarf objects. Our anal- 
ysis therefore demonstrates the importance of the low- 
luminosity galaxies, H i clouds and dwarf objects in the 
evolution of the galaxies. Interactions with dwarf galax- 
ies also may initiate star-formation events in normal spiral 



galaxies, such it occurs in the external arms of Mkn 1199, 
in Haro 15, surrounding Mkn 1087, o r in the i mpressive 
galaxy pair NGC 1512/1510 (Koribals ki fc Lope z-Sanchea 
2009) . Definetively, interactions between dwarf galaxies is 
one of the main triggering mechanism of the star- formation 
activity in starburst galaxies, but these dwarf objects are 
only detected when deep optical images and spectroscopy 
and complementary H i observations are obtained. 



12. Conclusions 

We have presented a comprehensive analysis of a sample of 
20 starburst galaxies that show the presence of a substan- 
tial population of very young massive stars, most of them 
classified as Wolf-Rayet galaxies. In this paper, the last of 
the series, we analyze the global properties of our galaxy 
sample using all multiwavelength data, that include X-ray, 
FUV, FIR, and radio (both Hi spectral line at 21 cm and 
1.4 GHz radio-continuum) results. Our main conclusions 
are the following: 

1. We compared the values of the Sfr derived from sev- 
eral indicators that consider fluxes at different wave- 
lengths. The results agree well within the experimental 
errors and with our Ha-based values, that were obtained 
after correcting for reddening and [Nil] contribution. 
However, we consider that t h e new Ha-based calibration 
provided by ICalzetti et al.l (|2007l) should b e preferred 
over t he well-known and extensively used iKennicufrg 
(|1998( ) calibration. Additionally, we checked that the 
FUV-h&sed Sfr very often shows similar results to 
those obtained using the emission of the ionized gas, 
providing a powerful tool to analyze independently the 
star-formation activity in both global and local scales. 

2. We checked that the Sfr/Lb ratio decreases with in- 
creasing metallicity. We derived empirical relationships 
between the [/-band, i3-band, and X-ray luminosities 
and the Sfr, that can be only used in starburst galax- 
ies and as a first estimation of the real Sfr value. 

3. All objects except one in our galaxy sample satisfy 
the FIR-radio correlation, indicating that they are pure 
star- forming systems. Only the galaxy ESO 566-8 lies 
away from the FZR-radio correlation because it seems 
to host some kind of nuclear activity. The non-thermal- 
to-thermal ratio seems to increases with increasing lu- 
minosity, suggesting that the cosmic-ray confinement is 
more efficient in massive galaxies than in dwarf objects. 

4. We provided empirical relationships between the ionized 
gas mass, neutral gas mass, dust mass, stellar mass, and 
dynamical mass with the i?-luminosity. Although all 
mass estimations increase with increasing luminosity, we 
find important deviations to the general trend in some 
objects, that seem to be consequence of peculiarities in 
these galaxies. The comparison between the dynamical 
mass (derived from the kinematics of the neutral gas) 
with the Keplerian mass (obtained from the kinematics 
of the ionized gas) and the stellar mass (from the H- 
band luminosity) provides further clues about systems 
in which the dynamics seem to be highly perturbed. We 
remark the importance of this study, because it is not 
common to find in the literature a comprehensive and 
detailed analysis of a sample of galaxies for which the to- 
tal (dynamical or stellar) mass, the reddening-corrected 
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luminosity in optical and NIR filters, and the To-based 
oxygen abundance, have been derived in a coherent way. 

5. We investigated some mass-metallicity relations and 
compared with pr evious results found in th e literature. 
As pointed out by iKewlev fc EllisonI ()2008[ ). the choice 
of the metallicity calibration has a strong effect in the 
derived M — Z relation. The tightness of the Mdyn — Z 
calibration indicates that that the dark matter con- 
tent also increases with metallicity. The scatter in the 
Mstars — Z and M^ar — Z relations are consequence of 
both the nature (dwarf galaxies, TDG candidates, merg- 
ers) and the star-formation histories experienced in each 
galaxy. 

6. We found that our sample galaxies agree well with 
the Sch midt-Kennicutt scaling l aw of star- formation de- 
rived bv lKennicutt et aP ([20071), that considers individ- 
ual star-forming regions within M 51. Some important 
deviation are found in NGC 5253 and Mkn 1199, that 
are very Hi-deficient, and in SBS 1319-1-579, where the 
star-formation activity seems to be supressed. 

7. The study of the mass-to-light ratios reinforces some of 
the results found in our analysis. We found that the 
neutral-gas-mass-to-luminosity ratio clearly decreases 
with increasing mass, as it seems to happen with 
the ionized-gas-mass-to-luminosity ratio. The ionizing- 
cluster-mass-to-luminosity ratio, however, seems to be 
constant with metallicity. The fact that we do not find 
any dwarf galaxy high Mm/ Lb ratio indicates that they 
have not experienced a lonely life. The analysis of the 
Mgas/ Mstars ratio suggests that this kind of galaxies 
have equal amount of neutral and stellar masses for 
metallicities 12-hlog(0/H)-8.2-8.3. The stellar-mass- 
to-luminosity ratio clearly increases with the B — R 
colour. 

8. We found that the reddening coefficient derived from 
the Balmer decrement clearly increases the the warm 
dust mass, indicating that the extinction is mainly in- 
ternal to the galaxy and not in the line-of-sight. We 
confirmed that the dust-to-gas ratio increases with the 
metallicity, and suggested that the low M dust /M gas ra- 
tios in dwarf low-metallicity galaxies is consequence of 
the large reserves of un-enriched neutral gas. However, 
the low M dust /M gas ratios observed at the outskirts of 
spiral galaaxies seem to be a result of the decreasing of 
the star-formation efficiency in these regions. 

9. The comparison of our data with the closed-box model 
clearly indicates that environment effects have played 
and important role in the evolution of the analyzed 
galaxies. The main effective yield we derived for our 
data agrees quite well with results found in the litera- 
ture, in particular with results found in other starburst 
or irregular dwarf galaxies. 

Considering all available data, we quantified how many 
galaxies are experiencing interaction or merger processes. 
We found that 17 up to 20 objects are clearly interacting 
or merging with low-luminosity dwarf objects or H i clouds, 
and all the remnant three galaxies (Mkn 5, SBS 1054-1-364, 
and SBS 1415-f437) show considerable divergences of some 
properties when comparing with similar objects. However, 
the interacting/merging features are only detected when 
deep optical spectroscopy and a detailed multi-wavelength 
analysis, remarking analysis of the kinematics and distri- 
bution of the neutral gas, are obtained. We therefore con- 



clude that interactions do play a fundamental role in the 
triggering mechanism of the strong star-formation activity 
observed in dwarf starburst galaxies. This observational re- 
sult completely agrees with the hierarchical model of galaxy 
formation that considers that large galactic structures were 
built up from the accretion of dwarf galaxies. 
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Appendix A: Final summary of individual galaxies 

In this Appendix we compile the main results found in 
our multiwavelength analysis of each individual system 
within our sample of Wolf-Rayet galaxies. In Papers I 
and II we described the optical/ A'liZ broad-band and 
Ha photometry, and the intermediate-resolution opti- 
cal spectroscopy analysis, respectively, of 16 up to 20 
galaxies studied in this work. The analysis of the 
additional 4 systems were presented in previous pa- 
pers: NG C 1741 (member AC with in the HCG 31 
group) in iLopez-Sanchez et al.1 (^2004a^■ Mkn 1087 an d 
their surrounding galaxies in lLopez-Sanchez et al.l (|2004bl) . 
the lum inous blue compact gal axy IRAS 08339 -1-6517 
in iLopez-Sanchez. Esteban fc Garcia-RoiasI ()2006[ ). and 
NGC 5253 in iLopez- Sanchez et al.l (|2007t) . Paper III com- 
piled the localization of the WR-rich star clusters within 
the galaxies and the analysis of their massive stellar pop- 
ulations. Paper IV compiles the global analysis of colours, 
and the physical properties and chemical abundances of the 
ionized gas. This paper, the last of the series, completes 
our analysis via a multi-wavelength analysis involving X- 
ray, UV, FIR, and radio data in both the 21-cm Hi line 
and the 1.4 GHz radio-continuum. 

— NGC 1741 hosts a very strong star-formation event, 
that is very probably a consequence of the merging 
of two spiral galaxies. NGC 1741 is the main member 
(AC) of the galaxy group HCG 31, and it is interacting 
with other galaxies in the group, including Mkn 1090 
(HCG 31 G). We detect both the blue and red WR 
bumps in its brightest region, as well as the nebular 
Hen A4686 hue. HCG 31 AC seems to have a slightly 
higher N/0 ratio. Some dwarf objects (members E, Fl, 
F2 and H) are tidal dw a rf gala xy (TDG) candidates. 
See lLopez-Sanchez et al.l (|2004al ) for details. 

— Mkn 1087 is a Luminous Compact Blue Galaxy 
(LCBG) in interaction with the nearby galaxy 
KPG 103a and with a dwarf surrounding galaxy (N 
companion). Deep optical images show long stellar tails 
connecting the main body of the galaxy with diffuse ob- 
jects, some of them hosting star-formation activity, and 
several TDG candidates. Although WR features were 
previous ly rep orted b y other a uthors , we do not detect 
any. See lLopez-Sanchez et al.l (|2004bD for details. 

— Haro 15 probably is a medium-size Sc spiral in inter- 
action with two nearby dwarf objects. Knot A shows a 
very high star-formation activity and WR features; its 
disturbed kinematics suggests that it is experiencing a 
minor merger with Haro 15. Knot B is an independent 
object because of its morphology, decoupled kinematics 
and chemical abundances. 

— Mkn 1199 is a system composed by a Sb-Sc spiral and 
a dwarf galaxy, both in clear interaction, as they may 
be at the first stages of a minor merger. The interaction 
has triggered the star formation activity in some areas 
of the main galaxy. Both the blue and red WR bumps 
are detected in the central region of Mkn 1199, that 
has solar metallicity. It seems that a substantial frac- 
tion of the H I gas has been expelled to the intergalactic 
medium because of its low Mm/ Lb ratio. 

— Mkn 5 is a BCDG with a strong star-forming burst lo- 
cated in the external part of the galaxy. The blue WR 
bump is detected in this starbursting region, that also 
possesses an important underlying old stellar compo- 



nent. We do not find any evidence of interaction, but 
the amount of H i gas of the galaxy is very low com- 
pared with that expected for a BCDG. Furthermore, 
its dynamical mass is higher than that expected for an 
dwarf galaxy with similar properties. Both results sug- 
gest that Mkn 5 has lost its neutral gas in some moment 
of its past and still has a disturbed H i kinematics. 
IRAS 08208+2816 is a Luminous Infrared Galaxy 
(LIRG) showing two long tails with a very high star- 
formation activity. The kinematics of the ionized gas 
clearly indicate merger features and the existence of 
two long tidal tails with TDG candidates. The chem- 
ical abundances of the brightest knots also seem to be 
different. We detect both the blue and red WR bumps 
in the central region, that possesses a high N/0 ratio. 
IRAS 08339+6517 is a luminous infrared galaxy 
(LIRG) and a luminous blue compact galaxy (LCBG) 
in clear interaction with a nearby dwarf galaxy. The 
majority of the H i gas of the system has been ex- 
pelled to the intergalactic medium because of this in- 
teraction (Cannon et al. 2004). Our deep optical im- 
ages reveal a faint stellar plume coincident with the 
H I tail, and a disturbed morphology in the outskirts 
of the galaxy. A particular bright knot may be a TDG 
candidate of the remnant of a previous minor merger. 
We detect weak WR features in its central burst and 
quantified the star formation history of the galaxy 
(iLop ez-Sanchez. Esteban fc Garcia-Roiasll2006l ). 
POX 4 is a morphology-disturbed low- metallicity 
BCDG showing strong star-forming bursts throughout 
all the galaxy. It seems to be in interaction with a 
nearby dwarf object that may have passed through the 
main body o f the galaxy, being the orig in of its ring-like 
morphology (jMendez fc Estebanlll997 ) and kinematics. 
However, this objects may also be a TDG candidate 
originated by the interaction with a n earby and diffuse 
H I cloud (|L6pez-Sanchez et al.ll20T0bl ). The He ii A4686 
emission line is clearly detected in its brightest region, 
as well as both the blue and the red WR bumps. 
UM 420: is a blue compact galaxy, but not a dwarf 
object, hosting intense star- formation activity. Besides 
it is located at 237 Mpc, we observe a central region 
and two kind of bright Ha tails pointing towards differ- 
ent directions. Its kinematics is also perturbed. It has a 
very low metallicity for an object with its absolute op- 
tical/ A^/i? luminosities, suggesting that it is a merging 
of two independent galaxies. We detect the He II A4686 
emission line but no the blue WR bump in its bright- 
est region. We found a probable N/0 enrichment in the 
central region. Its colors and properties are somewhat 
contaminated by the spiral disk of the foreground galaxy 
UGC 1809, located at 97 Mpc. 

SBS 0926+606: is a galaxy pair with high star forma- 
tion activity. Member A is a BCDG that shows a double 
nucleus; both its morphology and kinematics strongly 
suggests that it is a galaxy merger. We do not detect 
WR features in this galaxy but only the Hell A4686 
emission line. On the other hand, member B (another 
BCDG) hosts less star-formation activity, but it also 
shows hints of interactions, remarkably a long diffuse op- 
tical tail that shows a TDG candidate. SBS 0926+606 B 
has a huge emission in UV; the Sfr derived from the 
FUV luminosity is more than one order of magnitude 
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higher than the Ha-based Sfr. The system still hosts a 
huge amount of neutral gas. 

SBS 0948+532: is a very compact and blue object that 
hosts a very high star formation activity. Its A/hii/^b 
ratio is very high in comparison with objects with sim- 
ilar properties. Although usually classified as BCDGs, 
its total i3-luminosity indicates that it is not a dwarf 
object. We detect a faint optical tail mainly composed 
by old stars and with a slightly disturbed kinematics. 
We observe the nebular and broad He ii A4686 lines. 
SBS 1054+365: is a very nearby BCDG showing sev- 
eral star-forming regions embedded in a elliptical enve- 
lope composed by old stars. The kinematics of the ion- 
ized gas seems to be slightly disturbed. The main star- 
bursting region shows the nebular and broad He ii A4686 
lines. Although we do not detect any clear interaction 
feature, its dynamical mass is too high in comparison 
with that observed in similar objects, and its metallic- 
ity is too high for a dwarf object. Further studies are 
needed to clarify its nature. 

SBS 1211+540: is a very low-metallicity BCDG. It is 
composed by two bright Ha regions surrounded by a 
relatively old stellar component. This BCDG seems to 
show a higher metallicity than expected for a dwarf ob- 
ject with its same properties. The detection of two faint 
optical tails and its disturbed kinematics suggest that 
this galaxy is experiencing its first stages of a merger 
process. Although reported previously, we do not detect 
any WR features. 

SBS 1319+579: is a cometary-like BCDG showing two 
chains of intense star- forming regions over an underlying 
low-luminosity component dominated by old stars. We 
detect a very faint blue WR feature in the brightest 
knot. The analysis of the kinematics of the ionized gas 
strongly suggests that it is composed by two objects 
in interaction, that it is happening edge-on. Although 
there is plenty of neutral gas, the star formation is not 
very efficient, showing very low M-hh/Lb and Mi,/ Lb 
ratios in comparison with similar objects. Furthermore, 
it does not satisfy the Schmidt-Kennicutt law of star 
formation and the H i dynamics seem to be perturbed. 
We consider that the neutral gas has been expelled from 
the galaxy, but interferometric observations are needed 
to probe it. 

SBS 1415+437: is a very low-metaUicity BCDG that 
hosts a very strong star-forming region in which the neb- 
ular Hell A4686 emission line is observed. It possesses 
an important old stellar population underlying the star- 
burst. We do not detect any optical nearby companions 
and it does not show any evidence of interactions. 
Ill Zw 107: is a BCDG showing two strong star- 
forming bursts embedded in an irregular envelope. A 
diffuse prominent tail is detected in this object. The 
broad Heii A4686 line is found in the brightest knot, 
that shows a slightly higher N/0 ratio. The neutral gas 
may have been expelled and/or dispersed. This galaxy 
is likely composed by two dwarf objects in process of 
interaction or merging. 

Tol 9: is a BCG that belongs to the Klemola 13 galaxy 
group. It is a elliptical-shaped galaxy with intense neb- 
ular emission and chemically evolved. We have detected 
morphological and kinematical pattern that suggest in- 
teraction features. Our deep Ha image reveals an ex- 
tended filamentary structure with two main features 



that are located almost perpendicular to the main opti- 
cal axis of the galaxy. The probable origin of this struc- 
ture is a galactic wind. We detect both the blue and red 
WR bumps in the central region. The H i morphology 
and kinematics is quite intriguing, because this galaxy 
and two surrounding dwarf objects are embedded in the 
same Hi cloud (LS08b,LS+10b). 

Tol 1457-262: is a system composed by two bright ob- 
jects and two dwarf galaxies, all showing nebular emis- 
sion. We detect the nebular Hen line in the brightest 
knots of the main object. The regions within this sys- 
tem show chemical differences and peculiar kinematics. 
The neutral gas content seems to be very high, and its 
dynamics highly perturbed, although detailed H i map 
should be required to quantify this. We consider that 
this system is a galaxy group in which its members are 
in interaction. 

Arp 252: is a galaxy pair composed by two spiral galax- 
ies, ESO 566-8 (A) and ESO 566-7 (B), in the first stages 
of a major merger. This object shows two long tails 
mainly composed by old stars but hosting some star- 
forming regions and TDG candidates. ESO 566-8 shows 
the broad and nebular He ii A4686 emission line and 
the red WR bump. Its N/0 is quite high for a galaxy 
with its oxygen abundance. ESO 566-8 has a strong 
star- formation and may host some kind of nuclear activ- 
ity, because the FIR/ia.dio relation is not satisfied in it. 
Although it was previously observed by other authors, 
we do not observe any WR feature in ESO 566-7. 
NGC 5253: is a very nearby BCDG showing many 
peculiarities with respect to objects of similar charac- 
teristics. We detect clear broad WR features in the cen- 
tral regions, indicating the presence of both WNL and 
WCE stars. We confirmed the presence of a localized N 
enrichment in certain zones of the center of the galaxy 
and suggested a possible slight He overabundance in the 
same areas. We demonstrated that the enrichment pat- 
tern agrees with that expected for the pollution by the 
ejecta of WR stars. The amount of enriched material 
needed to produce the observed overabundance is con- 
sistent with the mass lost by the number of WR stars 
estimated in the starbursts (see Lopez-Sanchez et al. 
2007 for details). Although the kinematics of the ion- 
ized gas is somewhat peculiar, stellar kinematics seem 
to be consequence of rotation. Our optical study has not 
reveal any disturbed feature of a recent interaction pro- 
cess. However, its H i morphology is disturbed and its 
kinematics is quite intriguing, because it does not show 
any sign of regular rotation. The origin of this anomaly 
is most likely the disruption/accretion of a dwarf gas- 
rich companion or the i nteraction with another galaxy 
in the M 83 subgroup (iKobulnickv fc Skillmanri2008l: 
iLopez-Sanchez et al.l l2008al l2010aD . Furthermore, its 
Mm/ Lb and M^ust/ Lb ratios are very low and it does 
not satisfy the Schmidt-Kennicutt law of star formation. 
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